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ABSTRACT 

We present Far Ultraviolet Spectroscopic Explorer observations of the interstellar gas toward 30 white 
dwarf and 1 subdwarf (SdO) stars. These sightlincs probe the Local Bubble (LB) and the local interstellar 
medium (LISM) near the LB. We systematically measure the column densities for the following species: 
C II, C II*, C III, N I, N II, O I, Ar I, Si II, P II, Fe II, Fe III, and H 2 . Our survey detected only 
diffuse H 2 molecular clouds (/h 2 < 10 -4 ) along six sightlines. There is no evidence from this study 
that H 2 exists well inside the perimeter of the LB. The kinematical temperature for H 2 is less than the 
usual temperature observed in the local interstellar clouds, implying different gas phases in the LISM. 
The relative abundance ratios of Si II, P II, and Fe II give insight about the dust content. These ratios 
vary, but are similar to the depletion patterns observed in warm and halo diffuse clouds in more distant 
sightlincs in the Galaxy. The N I/O I and Ar I/O I ratios are significantly subsolar within the LB. Outside 
the LB a larger scatter is observed from subsolar to solar. Because Ar and N are only weakly depleted 
into dust grains if at all, the deficiencies of their neutral forms are likely due to photoionization. The 
evidence for significant ionization of N (and hence Ar) is strengthened by the detection and measurement 
of N II, which is a dominant ion for this element toward many sightlines. C III appears to be ubiquitous 
in the LISM toward our sightlines, but C II remains the dominant ionization stage of C. The limits on 
Fe III/Fe II imply that Fe II is the dominant ion. These observations imply that photoionization is the 
main ionization mechanism in the LISM and do not support the existence of a highly ionized condition in 
the past. In view of the variations observed in the different atomic and ionic ratios, the photoionization 
conditions vary significantly in the LB and the LISM. The cooling rate in the LISM, l c (in ergs" 1 (H I 
atom) -1 ), derived from the emission of the C II 157.7 /zm line has a mean value of log l c — ~ 25.59±0.32 
dex, very similar to previous determinations. 

Subject headings: ISM: clouds - ISM: structure - ISM: atoms - ISM: abundances - ultraviolet: ISM - 
white dwarfs 



1. INTRODUCTION 

Because of its proximity, the local interstellar medium 
(LISM) provides an unique opportunity to study in de- 
tail the physics of the warm (partially ionized) interstellar 
gas. Such gas is a major component of the interstellar 
material in the Galaxy and in other galaxies. The overall 
characteristics (temperatures, densities, kinematics) of the 
LISM were determined with ground-based telescopes and 
the Copernicus, International Ultraviolet Explorer (IUE), 
the Hubble Space Telescope (HST) satellites (e.g., McClin- 
tock et al. 1975; Bruhweiler & Kondo 1982; Frisch & York 
1983; Linsky et al. 1993; Lallement et al. 1995; Redfield et 
al. 2002). Measurements of interstellar gas after it has en- 
tered the heliosphere has also provided information on the 
interstellar medium very close to the sun (see, e.g., studies 
of the local interstellar cloud by Lallement 1998 and by 
Slavin & Frisch 2002 and references therein). 

From the ground-based observations, via mainly the 
Ca II K and Na I D absorption lines, very accurate ab- 
solute wavelength and very high spectral resolution pro- 
vide precise information on the local gas dynamics and the 



complex velocity structure of nearby gas (e.g., Vallerga et 
al. 1993; Ferlet 1999, and references therein). UV reso- 
nance absorption lines are much stronger than the avail- 
able optical absorption lines (which are not in the dom- 
inant ionization stages), and are therefore more sensitive 
probes of the warm, low density LISM. Major advances 
in determining the physical conditions have been carried- 
out with UV observations, especially in the last decade, 
with the Goddard High Resolution Spectrograph (GHRS) 
and the Space Telescope Imaging Spectrograph (STIS) on- 
board HST (e.g., Linsky et al. 1993; Lallement et al. 1995; 
Lallement 1996; Linsky 1996). 

With high spectral resolution ultraviolet observations, 
the physical properties of the cloud surrounding our so- 
lar system (the local interstellar cloud, LIC) can be esti- 
mated. The LIC is not very dense («h ~ 0.1 cm -3 ), but it 
is warm (T ~ 7000 K) and partially ionized (e.g., Frisch & 
York 1983; Lallement 1996; Linsky 1996). Near the LIC, 
there are similar clouds with different velocities (see, Lin- 
sky et al. 2000). All these clouds constitute collectively 
the LISM, and the clouds within roughly 100 pc are in a 
region called the Local Bubble (LB; the cavity has a radius 
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of about 100 pc in almost every direction, except toward 
one very low column density direction (/? CMa) where it 
extends to ^200 pc; Sfeir et al. 1999). Na I absorption line 
studies by Sfeir et al. (1999) show that the boundary of the 
LB is delineated by a sharp gradient in the neutral gas col- 
umn density with increasing radius. The distance of this 
boundary deduced from X-ray data is generally smaller 
by about 30% (see Snowden et al. 1998). In this study, 
we use the Far Ultraviolet Spectroscopic Explorer (FUSE) 
satellite to survey sightlines with distances less than 200 
pc. For the purposes of this study, this distance will define 
the extent of the LISM. Even though the distances are not 
great in the LISM, the sightlines are expected usually to 
have complex kinematical structure (e.g., Lallcment et al. 
1995). 

EUVE observations also have improved our knowledge 
of the ionization structure of the LISM (see the review by 
Bowyer, Drake, & Vennes 2000, and references therein). 
The EUVE data show that the neutral helium to neutral 
hydrogen ratio is about 0.07, which is somewhat less than 
the cosmic abundance of He to H of 0.1 (Dupuis et al. 
1995). Helium is therefore more ionized than H. Also, He 
shows less variability in its degree of ionization from one 
region to the next (Wolff, Koester, & Lallement 1999). 
One possible explanation for the high fractional ionization 
of He in the LISM is that there has not been enough time 
for the ionized He to recombinc to an equilibrium concen- 
tration from a more highly ionized condition in the past, 
perhaps from the influence of radiation from a supernova 
or its Shockwave less than few 10 6 years ago (Frisch & 
Slavin 1996; Lyu & Bruhweiler 1996). A different alter- 
native is that the LISM is currently exposed to a strong, 
steady flux of photons with E > 24.6 eV (e.g., a conductive 
interface between the LISM and a surrounding medium at 
T ~ 10 6 K; Slavin & Frisch 2002, or by recombination 
radiation from highly ionized but cool gases; Breitschw- 
erdt & Schmutzler 1994). A way to distinguish between 
these possibilities is to study the lightly depleted elements 
Ar, N, and O (Sofia & Jenkins 1998). In particular, the 
ionization of Ar (and to a lesser extent of N) can help us 
to differentiate between photoionization equilibrium and 
non-equilibrium cooling models. 

It is possible to increase our knowledge of the ioniza- 
tion structure of the LISM significantly by using data from 
FUSE (Jenkins et al. 2000; Lehner et al. 2002; Moos et al. 
2002; Wood et al. 2002b). For instance, Jenkins et al. 
(2000) found a deficiency of Ar I toward 4 white dwarf 
(WD) stars, thus favoring photoionization equilibrium. In 
the present work, we increase the sample to 31 stars that 
are located within 200 pc. WDs are particularly good 
background objects for studying the LISM because they 
are nearby, they have relatively simple stellar continua, 
and they are often nearly featureless. FUSE gives access 
to the wavelength interval between 905 A and 1187 A, a re- 
gion containing features of some important neutral species 
(N I, O I, and Ar I), many ionized species (C II, C III, N II, 
P II, Fe II, and Fe III), and H 2 . Surveys of O VI interstellar 
absorption and D I/O I ratio will be discussed by Oegerle 
et al. (2003) and Hebrard et al. (2003), respectively. Al- 
though FUSE has a moderate resolution of R ~ 20, 000 
and, hence, only the properties of the dominant cloud or 
the "average" properties of the clouds can be derived, the 



large number of resonance lines for atoms in several ioniza- 
tion stages as well as molecular lines provides the unique 
opportunity to study in detail the ionization structure, 
dust, and molecules in the LISM. 

2. OBSERVATIONS 
2.1. The Sample 

In Table 1, we list the positions, spectral types, dis- 
tances, and heliocentric velocities for the sample of 31 stars 
observed with FUSE. The distribution of the targets on the 
sky is shown in Galactic coordinates in Figure 1. There is 
a good coverage of the sky in almost each quadrant with 
the exception of the region spanned by I = 180° — 360° 
and b > 0°. Most of the WDs are from the FUSE PI 
Team programs (P104 and P204) used to study D/H and 
O VI in the LISM. The primary criterion for selecting the 
stars in Table 1 was that we would be able to measure the 
species of interest for this work (see § 3) with a sufficient 
degree of confidence. 

Recently, in the course of the FUSE observatory pro- 
gram "Survey of White Dwarfs from the McCook and Sion 
Survey" (Z903), several other WDs were observed. The 
quality of these data are not good enough for this study 
because strong airglow emission lines are present at all 
wavelengths and the signal-to-noise levels are too low (the 
flux of these WDs in this program are systematically lower 
than in the PI team programs). We also note that these 
targets do not fill in the deficiency of targets in the quad- 
rant at I = 180° - 360°, b > 0°. 

2.2. Instrument and data reduction 

The FUSE instrument consists of four channels: two 
optimized for the short wavelengths (SiC 1 and SiC 2; 905- 
1100 A) and two optimized for longer wavelengths (LiF 1 
and LiF 2; 1000-1187 A). There is, however, overlap be- 
tween the different channels, and, generally, a transition 
appears in at least two different channels. More complete 
descriptions of the design and performance of the FUSE 
spectrograph are given by Moos et al. (2000) and Sahnow 
et al. (2000). To maintain optimal spectral resolution 
the individual channels were not co-added together. Also, 
measurements on the independent spectra could be com- 
pared with each other to obtain information about the 
effects of detector fixed pattern noise. Detailed informa- 
tion for each sightlinc, including whether the low resolu- 
tion (LWRS) apertures or the medium resolution (MDRS) 
apertures were used, can be found in Table 1. 

Standard processing with the current version of the cal- 
ibration pipeline software (version 2.0.5 and higher) was 
used to extract and calibrate the spectra. The software 
screened the data for valid photon events, removed burst 
events, corrected for geometrical distortions, spectral mo- 
tions, satellite orbital motions, and detector background 
noise, and finally applied flux and wavelength calibrations. 
The extracted spectra associated with the separate expo- 
sures were aligned by cross-correlating the positions of ab- 
sorption lines, and then co-added. In some cases, the lack 
of strong absorption lines or the contamination of airglow 
lines prevented us from determining the relative displace- 
ments, so we simply coadded the spectra with no shifts. 
Generally, the different exposures were obtained succes- 
sively, so no large shifts were expected and no significant 
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loss of spectral resolution was observed. The co-added 
spectra were finally rebinned by 4 pixels (~ 27 mA) since 
the extracted data are oversampled. 

Figure 2 presents two examples of FUSE calibrated spec- 
tra of WDs. As can be seen in the figure, the stellar con- 
tinua are relatively simple and can be fitted with low-order 
(< 3) Legendre polynomials. The primary error sources 
are statistical noise, fixed-pattern noise, and possible sat- 
uration effects in the lines. The fixed-pattern noise is re- 
duced by comparing the same spectral feature in different 
channels, while saturation can be examined when there 
are several transitions of the same species with different 
/A values (§ 3). 

3. SPECTRAL FEATURES 
3.1. General Considerations 

The primary objective of this work is to measure the col- 
umn densities of the neutral atoms, weakly ionized atoms 
and molecular hydrogen toward sightlincs within ^200 pc 
in order to probe the physical conditions within the LISM. 
For such study FUSE provides access to several important 
resonance lines, and we have measured systematically the 
column densities, N, (or limits on the column densities) 
for the following species: C II, C II*, C III, O I, N I, N II, 
Si II, P II, Ar I, Fe II, and Fe III. The column densities 
of H 2 were determined by using the Lyman and Werner 
bands. We note that P III and N III also are present in the 
FUSE bandpass and are potentially interesting ions. How- 
ever, the P III transition is too weak to be detected and 
its 3a upper limit does not constrain the ratio P II/P III 
in a meaningful way. N III is blended with the saturated 
Si II line and possibly with photospheric N III making any 
measurements of this line highly uncertain. However, up- 
per limits obtained by assuming that N III is not blended 
and lies on the linear part of the curve of growth (COG) 
show that N III is not the dominant ion stage of N in the 
LISM. Corrections for saturation would increase the up- 
per limit, but this effect is likely less than the contribution 
from the Si II blend. Therefore, nitrogen is primarily N I 
and N II in the LISM. 

In Table 2, we summarize the atomic/ionic transitions 
used in this work. To measure the column densities (see 
below), we adopted wavelengths and oscillator strengths 
from the Morton (2000, private communication) atomic 
data compilation. This compilation is similar to the Mor- 
ton (1991) compilation with only a few minor updates to 
the atomic parameters for lines of interest in this study. 
For the Fe II lines, the new oscillator strengths derived 
by Howk et al. (2000) were adopted. We note, however, 
that there is a new laboratory measurement of the /-value 
for Fe II A1144.9 (/ = 0.083 ± 0.006, Wiese, Bonvallet, 
& Lawler 2002), which is about 30% lower than the ex- 
perimental value of Howk et al. (2000). However, Wiese, 
Bonvallet, & Lawler (2002) note that additional study is 
necessary for the other transitions available in the FUSE 
bandpass, rather than a simple rescaling, as they are fur- 
ther down on the COG than the Fe II Al 144.9 line. As we 
use all the transitions available in the FUSE spectrum, we 
employ the /-values of Howk et al. (2000). 

For molecular hydrogen, we use the wavelengths from 
Abgrall et al. (1993a, b). The H 2 /-values were calculated 
from the emission probabilities given by those authors. 



3.2. Photospheric Lines 

To varying degrees, WD stars have sharp metal photo- 
spheric lines that can mimic interstellar lines, and there- 
fore it is important to identify and disentangle them. Some 
WDs can be pure hydrogen (no photospheric metal lines) 
and some have only small amounts of metals in their pho- 
tosphere. Their metallicity is generally very sub-solar (less 
than a few 10~ 3 solar metallicity, Zq), and only a few 
metal-rich WDs are as high as 0.1 Z G (e.g., WD 2211-495, 
WD 0621— 376). As a consequence, WDs have relatively 
smooth continua with only a few stellar lines. Figure 2 
presents spectra of two WDs, one pure hydrogen and the 
other with some metals, with the interstellar lines flagged. 
Table 1 summarizes the metallicity of the WDs. 

The lines that can be affected by photospheric lines are 
C III A977, Ar I A1066, P II A1052, and Fe III A1122. C III 
can be present in WDs with 20 x 10 3 < T cff < 50 x 10 3 K, 
and peak at ~ 30 x 10 3 K (see Chayer, Fontaine, & We- 
semael 1995). Thus, when the star is not pure hydrogen, 
some contamination is expected for C III. Generally the 
stellar radial velocity is not different enough from the in- 
terstellar component to be able to distinguish between the 
two components. Photospheric Fe III can be also present 
for WDs with T off < 40 x 10 3 K. Moreover, Fe III A1122.53 
is blended with photospheric Si IV A1122.49. For this 
species, we have a direct diagnostic for knowing if there 
is photospheric Si IV A1122 or Fe III A1122 by inspect- 
ing the lines Si IV A1128 and Fe III A1128 (see Figure 
2 in Chayer et al. 2000). P II and Ar I are not found 
in the photosphere of the WDs, but can be blended with 
Fe VI, Ni VI A1152 and Si IV A1066, respectively. Fe VI 
(and Ni VI) are present in significant quantities only at 
T cff > 60 x 10 3 K (Chayer, Fontaine, & Wesemael 1995). 
Therefore for most of our sample, these blends with Fe VI 
and Ni VI are not much of a concern. For Si IV A1066 we 
look at longer wavelengths (e.g. Si IV A1128) for possible 
presence of this line. We finally note that some of the N I, 
O I, Fe II, and H 2 transitions may be blended with some 
stellar lines, but the (generally) large numbers of transi- 
tions for these species allow us to detect if a specific line 
is contaminated. 

3.3. Circumstellar Lines 

Our survey used hot WDs, typically with tempera- 
tures ranging from 30,000 to 70,000 K (the subdwarf SdO 
BD+28°4211 has a temperature of 82,000 K; Sonneborn et 
al. 2002). Such hot WDs are emitting a large number of 
ionizing ultraviolet photons, creating Stromgren spheres 
with radii from 8 upto 30 pc depending on the electron 
density and the temperature of the WD (Dupree & Ray- 
mond 1983; Tat & Terzian 1999). Bannister et al. (2003) 
recently found direct evidences of such circumstellar ab- 
sorption features (in form of high ions such as C IV) in 
the spectra of 4 WDs out of 11 present in our survey. 
Therefore, our study of the ionization conditions in the 
LISM could a priori be systematically biased because of 
the use of hot WDs. However, to probe these conditions, 
we employ mainly diagnostics from neutral species, such 
as N I, O I, and Ar I, species not present in the completly 
ionized Stromgren spheres. We have found no correlation 
between the ratios of Ar I/O I and N I/O I with the tem- 
perature of the WDs. While it is clear that the hot WDs 
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participate to the overall ionization of the LISM (Vallerga 
1998; Slavin & Frisch 2002), we have no evidence that 
stronger EUV sources give different results compared to 
weaker EUV sources in our column density ratios. 

3.4. Geocoronal Emission Lines 

Contamination of the spectra by geocoronal emission 
can be significant for the LWRS apertures of the FUSE 
spectrograph. Such emission lines are very strong at H I 
Ly [3, but decline rapidly in strength for the higher Ly- 
man lines. During the sunlit part of the orbit, airglow 
emission lines are also present for many of the same lines 
of N I, O I, and Ar I used to measure the absorption by 
interstellar species. The importance of airglow emissions 
depends on many factors, including the fraction of sunlit 
time, the limb angle, the flux of the target, the aperture 
used, and the length of individual exposures. For the stars 
considered in this study, the observational strategy was to 
reduce the contribution from airglow lines: (1) The tar- 
gets were observed in a manner that maximized the night 
time fraction, ranging typically from 30% to 100% (only 
WD 1234+481 was observed with basically no night time). 
(2) In order to remove small spectral shifts during an ob- 
servation caused by thermally induced motion of the op- 
tics, an exposure is broken up into a large number of sub- 
exposures. Correlation techniques using uncontaminated 
interstellar lines are used to remove the shifts between the 
spectra and they are then added. To the extent the shifts 
are due to telescope mirror motions rather than grating 
motions, these shifts will smear out the airglow contribu- 
tion. In practice, impact of the airglow contamination is 
reduced by this technique. (3) The objects studied are 
relatively bright with fluxes from a few 10~ 12 to 10~ n 
ergcm~ 2 s _1 A -1 . (In contrast the stars from the Z903 
program have lower fluxes, long individual exposure times, 
and are observed largely in day time.) However, even un- 
der the conditions mentioned above, the airglow could still 
be significant. Therefore, we took several steps discussed 
below to check for possible contamination. 

The widths of the airglow lines are determined by the 
size of the spectrograph apertures. For the LWRS aper- 
ture, the width of an airglow is typically 0.35 A, while it is 
typically 0.045 A for the MDRS aperture. Figure 2 illus- 
trates the differences in strength and width of the airglow 
lines between the MDRS aperture (WD 1800+685, night 
fraction of 0.85) and LWRS aperture (WD 0131-163, 
night fraction of 0.95). Hence with the LWRS aperture, 
the primary concern for measuring the interstellar lines 
considered here is the continuum placement, but for the 
MDRS aperture, the width is small enough that the emis- 
sion could fill the interstellar lines. Several of our targets 
were observed with both the LWRS and MDRS apertures. 
The column densities generally agree within the errors, 
and in particular there was no systematic smaller column 
densities derived with the MDRS compared to those de- 
rived with the LWRS. 

Feldman ct al. (2001) have reported the brightnesses 
of the terrestrial day airglow lines using the FUSE in- 
strument. We have examined airglow lines not blended 
with interstellar lines to estimate if our interstellar mea- 
surements could suffer from airglow contamination. At 
A > 1080 A, we use the O I 1 D- 1 D° A1152 line as an indi- 



cator for the presence of N I emission near 1134 A, since 
the lines have comparable strengths. This airglow line can 
cause a continuum placement problem for the measure- 
ment of P II A1152 as shown in Figure 2 in the spectrum 
of WD 0131-163. We were careful to use the N I A1134 
lines only when these lines were not contaminated by air- 
glow. At 1000 < A < 1080 A, we checked the strengths 
of several O I airglow lines at 1027, 1028 and 1042 A. 
These lines have strengths that indicate the strengths of 
the O I A1039 and Ar I AA1048, 1066 airglow lines. We 
do not find any cases where the contamination could be 
troublesome (in any case the ler errors quoted here should 
be larger than the effects of contamination). Finally, at 
A < 1000 A, where the most useful lines of N I and O I are 
located, the airglow emissions are negligible based on O I 
A990. Generally, we are confident that the airglow lines do 
not seriously compromise our interstellar measurements in 
ways we do not recognize. 

3.5. Interstellar Absorption Lines 
3.5.1. Species with Only One or Two Transitions 

The apparent optical depth (AOD) method (Savage & 
Sembach 1991) was used when only one or two lines of 
the same atomic species were available. The absorption 
profiles were converted into apparent optical depths per 
unit velocity, T a (v) — hn[I c / I \> s (v)}, where I c , I h s are the 
intensity without and with the absorption, respectively. 
r a (v) is related to the apparent column density per unit 
velocity, N a (v) (cm~ 2 (kms -1 ) -1 ) through the relation 
N a (v) = 3.768 x 10 U T a (v)/[f A(A)]. The integrated appar- 
ent column density is equivalent to the true integrated col- 
umn density in cases where the absorption is weak (r < 1) 
or the lines are resolved. This method was also used to 
place a lower limit on a saturated line. 

For C II and N II, only a strong transition is available, 
and therefore, except for a few cases, the column densities 
are lower limits. The P II transition is weak for all of the 
sightlines considered. C II* is also usually relatively weak 
for the sightlines considered, but can be slightly blended 
with the Lyman 5-0 R(l) line of H2. When hydrogen 
molecules are present along a sightline, the contribution 
of both lines can still be measured because the shift be- 
tween the two lines is large enough. 

Si II has two transitions. The weakest one is free of 
blending with other interstellar or stellar lines, but the 
strongest transition at 989.873 A is usually saturated and 
blended with N III of interstellar or stellar origin. 

For Ar I, two transitions are also available at 1048.22 
and 1066.66 A, with the latter being the weakest transi- 
tion of the two. We can therefore check for the effect of 
saturation, except when the WD is not pure hydrogen in 
which case Ar I A1066.66 is blended with the stellar Si IV 
A1066.66. 

3.5.2. Species with several (> 2) transitions: N I, Ol, 
Fe II, and H 2 

A COG analysis was undertaken when several lines were 
present (N I, O I, Fe II, and H2) by using the measured 
equivalent widths. The equivalent widths were measured 
by directly integrating the intensity, and their errors were 
derived following the method described by Sembach & Sav- 
age (1992). A single component Gaussian (Maxwellian) 
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curve of growth was constructed in which the Doppler pa- 
rameter b and the column density were varied to mini- 
mize the x 2 between the observed equivalent widths and 
a model curve of growth. For H2, first a COG was con- 
structed for the set of lines within each rotational level 
(J = 0, 3). The resulting column densities and b- values 
were then used as the starting point for a \ 2 minimization 
of a COG that simultaneously included all of the rotational 
levels with the column density for each level and a single 
b value as free parameters. 

Depending on the quality of the data and the amount 
of matter along a given sightline, we found that generally 
the column densities of N I, O I, and Fe II were well deter- 
mined with a COG because they have several transitions 
with different fX in the FUSE bandpass. We found that 
fixed-pattern noise was not a serious problem, and the ef- 
fects of saturation can be reasonably well understood. In 
Table 3, we compare the b- values of N I and O I obtained 
from the COG analysis, and within the errors, they are 
generally similar. In isolated cases where there was no 
good constraint on b, we used the 6-value from one species 
to constrain the column density of the other. Those cases 
are indicated by colons in Table 3. The different Fe II 
transitions generally lie on the linear part of the COG. 

Even though we considered seriously the saturation 
problem, with the FUSE spectral resolution, it still could 
be that narrow, unresolved components could be concealed 
by broader features. Jenkins (1986) showed, however, that 
column densities remain accurate to <~10% provided that 
the lines are not heavily saturated (central optical depth 
such as t < 2) and b and r do not have a bimodal 
distribution function. Three sightlincs (WD 0501+527, 
WD 1254+223, and WD 1314+293) of our sample were ob- 
tained with high spectral resolution using STIS. Redfield & 
Linsky (2003), but see also Vidal-Madjar et al. (1998) and 
Kruk et al. (2002), derived the column densities of several 
species considered in our survey for these sightlines, and 
their total column density measurements are in agreement 
within la with ours. This shows that FUSE measurements 
are reliable despite the moderate spectral resolution of the 
instrument. It also shows that the oscillator strengths of 
the species available in the FUSE wavelengths are also re- 
liable. 

3.5.3. 3(7 Upper Limits 

The 3a upper limits for the equivalent widths are defined 
as Wmin = 3a 5 X, where a is the inverse of the contin- 
uum S/N ratio per resolution clement SX = 0.1 A FWHM. 
The corresponding 3cr upper limits on the column density, 
presented in Table 4, are obtained from the correspond- 
ing equivalent width limits and the assumption of a linear 
curve of growth. 

3.6. Results from Other Sources 

Table 1 lists the sources of the column densities when 
they were derived from other investigations. For the sight- 
lincs that are reported in the ApJS special issue on the 
determination of D/H with FUSE (Hebrard et al. 2003; 
Kruk et al. 2002; Lehner et al. 2002; Lemoine et al. 2002; 
Sonneborn et al. 2002; Wood et al. 2002a), we comple- 
mented the measurements with additional measurements 
of the column densities for species considered here but 



not in those papers. In addition, for WD 1631+781 and 
WD 2004-605, the values reported here are a mean be- 
tween our methods (AOD and COG) and the profile fit- 
ting method reported by Hebrard et al. (2003). We note 
that our values agreed at 1 - 2a level. For WD 1211+332, 
WD 2247+583 WD 2309+105, and WD 2331-475, the col- 
umn densities are from Oliveira et al. (2003). For 
WD 0004+330, they are from Oliveira et al. (private com- 
munication, 2003). They were also obtained using several 
independent methods. In addition, we made independent 
measurements of the column species for WD 0004+330 and 
WD 2331-475; our values agreed at the 1 - 2a level. 

The column density results are summarized for each 
sightline in Table 4 for the atomic species and in Table 5 
for the molecular hydrogen. All the errors (except upper 
limits) arc la. 

Table 5 also lists the rotational excitation temperature 
of H2 for each sightline. A single straight line with a slope 
characterized by T ex fits all the observed levels J — 0, 3, 
except toward WD 0004+330 (see § 5). No rotational level 
higher than J = 3 was detected. 

4. GENERAL FINDINGS 

4.1. Sky-Distribution of the Ol Column Density 

O I is an excellent tracer of H I because its ionization 
potential and charge exchange reactions with hydrogen en- 
sure that the ionization of H I and O I are strongly cou- 
pled. Only a few sightlines considered in this work have 
an accurate H I column density measurement (see Moos et 
al. 2002). Measuring H I column densities is particularly 
complex for LISM sightlines because the H I Lyman series 
transitions in the FUSE wavelength range are generally 
on the flat part of the COG. In addition, the interstel- 
lar lines are narrow enough that blending with the stellar 
spectra can be a significant source of uncertainty. There- 
fore, we use O I as a proxy of H I, adopting the O I/H I 
ratio derived by Oliveira et al. (2003) (who considered one 
sightline plus the sightlines in Moos et al. 2002), O I/H I = 
(3.17 ± 0.19) x 10~ 4 , to estimate the amount of H I along 
a given sightline. This ratio is similar to values obtained 
for more distant sightlines (O I/H 1= (3.43 + 0.15) x 10~ 4 ; 
O I/H 1= (4.08 + 0.14) x 10~ 4 ; Meyer, Cardelli, & Sofia 
1997; Andre et al. 2003a, respectively), showing that O is 
usually a good proxy for H in different environments of the 
Galaxy. Also O is the most abundant species after H and 
He, making it a general good tracer of neutral gas over 
a large range of values of column density. The O I/H I 
ratio may vary in regions of elevated dust and density 
(Cartledge et al. 2001) or possibly in regions where the 
signature of stellar processed gas is still visible (Hoopes et 
al. 2003), but these conditions likely do not hold for the 
sightlines studied here. Therefore, we expect O I to be a 
reliable proxy for H I in this study. 

Figure 1 shows the distribution of the O I column den- 
sity (and thus H I) in the Galactic sky. The diameter 
of each circle is inversely proportional to the WD's dis- 
tance and the intensity of the shading of the symbol indi- 
cates the total column density along a given line of sight. 
Because the clouds are not resolved (see below), most of 
the sightlines have column densities between 15.5 and 16.5 
dex. Generally, the farthest sightlines result in the largest 
column densities. However, there are cases where, for ex- 
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ample at -90° < b < 0° and 300° < I < 360°, sightlines 
with very different stellar distances (stellar distances be- 
ing only an upper limit on the distances of the interstellar 
clouds) have a similar column density. At d < 100 pc, the 
North Galactic pole sightlines seem to have much lower 
column densities than those near the South Galactic pole 
(see also, Genova et al. 1990; Redfield et al. 2002; Welsh 
et al. 1999). The lowest column density is found near 
I ~ 190° (WD 0549+158) which is known for its relatively 
low density (Gry & Jenkins 2001; Redfield et al. 2002). 

4.2. Sightlines and the Local Bubble 

A typical sightline goes through the Local Interstellar 
Cloud (LIC), other warm partially ionized clouds, the Lo- 
cal Bubble (LB), and, depending on the distance of the 
star, might even go beyond the LB. With very high res- 
olution observations (R > 10 5 -10 6 ), the different clouds 
along a given sightline can be resolved as long as they do 
not occur at the same velocity. However, with FUSE, a 
resolution of only about 20 kms -1 is achieved. As a result, 
the column densities derived in this study reflect the total 
amount of matter along a given line of sight. 

Sfeir et al. (1999) have used the absorption line studies 
of Na I D, a tracer of cold neutral gas, to map the bound- 
ary of the LB. The LB consists of an interstellar cavity of 
low density neutral gas with radii between 65-250 pc in 
the galactic plane and ^100 pc in other directions which 
is surrounded by a denser neutral gas boundary (also re- 
ferred as the "wall"). The H I density corresponding to 
the wall is about 19.3 dex (20 mA equivalent width for 
Na I D). Therefore both the distance and total H I column 
density are used as a guide on whether or not a sightline 
extends beyond the edge of the LB. 

Based on these arguments, the clouds toward the follow- 
ing targets should be well inside the LB (log N(H I) < 19.0 
dex): WD 0050-332, WD 0455-282, WD 0501+527, 
WD 0549+158, WD 0621-376, WD 1254+223, 
WD 1314+293, WD 1634-573, WD 2111+498, 
WD 2152-548, WD 2211-495, and WD 2331-475. In 
contrast, the following sightlines have 19.0 < log7V(H I) < 
19.3 dex and therefore must be close to the bound- 
ary: WD 1211+332, WD 1615-154, WD 1636+351, 
WD 2004-605, and WD 2309+105. All the other sight- 
lines have logiV(H I) > 19.3 dex, though we note that 
some of them have values close to the approximate bound- 
ary. In short, some of our sightlines are well inside the 
LB, some are very close to the wall, and some are beyond. 
We note that the column densities in the LIC are of the 
order of 5 x 10 16 -2 x 10 18 cur 2 (Redfield & Linsky 2000). 
Therefore material in the LIC generally constitutes only a 
very small fraction of an observed total column density. 

We emphasize again that even though some of the sight- 
lines reside completely within the LB, it is likely that sev- 
eral clouds separated by less than 20 kms^ 1 are present 
along a given sightline. Therefore, with the FUSE obser- 
vations we only derive either the properties of the domi- 
nant (larger column density) cloud or a mixture of these 
clouds properties. This is the case, for example, for G191- 
B2B and GD 246 where 2 or 3 clouds are detected us- 
ing higher spectral resolution (Vidal-Madjar et al. 1998; 
Oliveira et al. 2003). This hidden component structure 
can subsequently complicate the physical interpretation of 



these spectra (see § 6 and 7). 

4.3. Distribution of the Column Densities for the 
Different Species 

In Figure 3, we present the column density distribution 
of O I, N I, Ar I, Si II, Fe II, and P II (based only on the 
measurements - not the limits - presented in Table 4) . Re- 
cently, Redfield et al. (2002) presented the column density 
distribution of individual clouds along sightlines within 100 
pc. A comparison with the distribution of Fe II column 
density in their work (their Figure 9) and this work clearly 
shows an increase in our sample of <~ 1 dex in the column 
density. Our targets generally lie at larger distances and 
thus our sightlines have intrinsically more gas along the 
line of sight. This difference is also due to the fact that 
our sightlines contain several clouds, which arc unresolved 
in the FUSE data. 

Figure 3 also gives some general information about the 
composition of the LISM. Not surprisingly, the highest col- 
umn densities are observed in the most abundant metals, 
O I and N I. These high column densities also indicate 
that the LB, known mostly for its hot gas that emits soft 
X-rays, nevertheless still has a substantial amount of neu- 
tral gas. Ar I appears to be deficient with respect to O I, 
suggesting that some of the gas is partially ionized (see 
§6). Si II and Fe II have similar solar abundances, but 
their column density distributions clearly show that Fe II 
is deficient with respect to Si II, suggesting that Fe is more 
depleted onto dust grains. P II is much less abundant than 
Fe II and Si II, but a comparison with the solar abundance 
(~ —2 dex with respect to Fe II and Si II) implies that Fe 
and Si are depleted onto dust grains, and hence dust in 
the form of silicate and iron grains must be present in the 
LISM (see § 7). 

5. MOLECULAR CLOUDS IN THE LISM 

5.1. Previous Studies 

Copernicus was the first satellite to probe the H 2 Ly- 
man and Werner lines in absorption (Spitzer et al. 1973). 
Spitzer, Cochran, & Hirshfeld (1974) presented H2 results 
toward 28 stars lying between 110 pc and 2 kpc (we cor- 
rected these distances with the parallaxes measured by 
Hipparcos). Seven of these stars are situated between 111 
and 166 pc and are presented on the Galactic sky in Fig- 
ure 1 with the symbol "C" . After Copernicus was decom- 
missioned, CO emission was the primary technique used 
to probe molecules in the LISM. Magnani, Blitz, & Mundy 
(1985) discovered a very large number of CO clouds at high 
galactic latitude \b\ > 20°, and based on statistical argu- 
ments estimated that these clouds (known as MBM clouds) 
are within 100 pc. Two of these clouds, known as MBM 12 
{(I, b) = (159°, -34°)) and MBM 16 ((I, b) = (172°, -38°)) 
clouds were considered until recently to be the nearest 
known molecular clouds to the Sun, at a distance of ~65 
pc and 60 < d < 95 pc, respectively (Hobbs, Blitz, & Mag- 
nani 1986; Hobbs et al. 1988). However, two recent inde- 
pendent estimates of the distance place the cloud MBM 12 
at 275-300 pc (Luhman 2001; Andersson et al. 2002), i.e. 
at a much larger distance than initially thought. These 
recent results could imply that the MBM clouds might 
be farther away. Dame et al. (1987) presented a CO sur- 
vey of the Galaxy at |6| < 30°, but the clouds studied 
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were estimated to be farther away, d > 140 pc. Trap- 
ero et al. (1996) also found two other molecular cloudlets 
in 12 CO emission in the Galactic plane within 120 pc, 
(/,&) - (149°, -5°). 

Direct evidence of H 2 in the LISM has been made pos- 
sible again with the recent launch of FUSE. Note that be- 
cause H 2 is observed in absorption, the distance to the star 
directly places a firm upper limit on the distance of any 
molecular cloud. Recently, Gry et al. (2002) used FUSE 
data to study H2 in three lines of sight passing through the 
Chamaeleon complex ((I, b) <~ (300°, —15°)), which has an 
estimated distance of 150 pc. 



5.2. Distribution of H2 in the LISM 

Although the FUSE instrument is sensitive to H 2 
columns down to low values, ~ 10 14 cm~ 2 , there appears 
to be detectable amounts of H2 along six sightlines out of 
31 in our sample (Table 5 summarizes the results). The 
sightlines have distances between 94 < d < 159 pc and lie 
at 60° < I < 120° and -30° < b < 40° (see Figure 1). We 
also examined the lower signal-to-noise data from the Z903 
program (see § 2.1). The only two stars in the program 
(at the time of this work) that show H2 in their spec- 
tra are WD 0421+740 ((I, b) = (136.13°, +16.81°), d = 257 
pc) and WD 1725+586 ((l,b) = (87.22°, +33.92°), d = 
114 ± 23). 7 There is, however, no H 2 detection toward 3 
other WDs lying in this region of the sky (WD 1820+580, 
WD 1845+683 [d = 125 pc], WD 1943+500 [d = 98 
pc]). We also note that H 2 was detected toward Feige 
110 at even higher latitude ((l,b) = (74.09°, -59.07°), 
d = 179+ If pc, Friedman et al. 2002). These diffuse H 2 
clouds roughly lie in the same direction as that of some 
of the CO detections. Finally, another star which was 
not used in our study of atomic and ionic column den- 
sities, WD 0441+467 ((l,b) = (158.49°, +0.47°), d = 110 
pc), has very strong H 2 transitions with a larger number 
of rotational levels compared to the other sightlines. This 
cloud could be related to the two CO clouds observed by 
Trapero et al. (1996, see below). 

Most or all of the H 2 detected in this study may be 
close to the boundary of the LB or outside it. Of the six 
sightlines in Table 5, one, WD 1636+351, has an O I col- 
umn indicating it is close to the inside of the LB wall, 
log(7V(H I)) = 19.2 dex. The others are 0.3 to 1 dex 
higher. It is also possible to compare detections and 
non-detections along neighboring sightlines. Some of the 
Copernicus sightlines (Spitzer, Cochran, & Hirshfeld 1974) 
that show H 2 are marked with a "C" in Figure 1. Four of 
them are close to WD 1615—154 where no H 2 was detected. 
As shown in Figure 1, WD 1800+685 (d = 159 pc), which 
is listed in Table 5, is near WD 1631+781 (d = 67 pc), 
which shows no detectable H 2 . Likewise, WD 2011+395 
(d = 141 pc) and WD 2247+583 (d = 122 pc), also 
listed in Table 5, are close to WD 2111+498 (d = 50 pc), 
which shows no H 2 . Although small-scale structure in the 
medium could cause these differences, it is likely that the 
molecular clouds reported here reside at distances between 
~60 pc and ~140 pc. Using the Na I D2 contour map 



of Sfeir et al. (1999) (their Figure 3), toward the general 
direction of WD 1615-154 (I w 350° - 0°) the edge of 
the LB is about 90 pc, while it is about 65 pc toward 
the general direction of WD 1631+781 and WD 2111+498 
[I « 90° — 110°). We are forced to conclude on the basis 
of both the O I column densities and distance considera- 
tions, that aside from WD 1636+351, which appears to be 
just at the wall, there is no evidence in the present data 
set for molecular clouds within the LB. We note, however, 
that Sfeir et al. (1999) found dense, cold clouds toward 
this direction at ^85 pc. 

5.3. Properties of the Observed H 2 

The total H 2 column densities are modest ranging from 
1.7 x 10 14 cm- 2 toward WD 2011+395 to 2.4 x 10 15 cm" 2 
toward WD 2148+286. Toward the other sightlines in our 
sample, the H 2 amount is less than ~ 10 14 cm~ 2 . These 
column densities are generally smaller than the ones de- 
rived by Spitzer, Cochran, & Hirshfeld (1974) (2 sight- 
lines have log N(H 2 ) between 14.5 and 15.3 dex, the others 
marked in Figure 1 have logA^(H 2 ) > 19 dex). 

The molecular fraction is defined as /h 2 = 
2iV(H 2 )/LV(H I) + 2N (H 2 )]. To obtain the values for 
the H I columns, we use O I as a proxy, except toward 
WD 2148+286 which has a good H I determination (Son- 
neborn et al. 2002). The molecular fractions obtained 
range between 5.0 x 10~ 6 and 1.4 x 10~ 4 and are com- 
parable to values obtained for other Galactic sightlines 
with log[7V(H I) + 27V (H 2 )] < 20 dex (Savage et al. 1977). 
We note that these sightlines do not have higher depletion 
onto dust grains compared to the other sightlines (sec § 7) . 

In contrast, Trapero et al. (1996) found Af(H 2 ) ~ 
N(R I) ~ few x 10 20 cm" 2 toward Perseus ((l,b) ~ 
(149°, -5°), d ~ 120 pc) from 12 CO measurements. The 
FUSE data, for WD 0441+467, discussed in § 5.2, indicates 
a large amount of H 2 probably at a level of a few times 
10 19 -10 20 cm~ 2 . 8 The stars observed by Copernicus are 
also at similar distances and also have H 2 column densities 
between 10 19 -10 20 cm" 2 (Spitzer, Cochran, & Hirshfeld 
1974). This indicates that the LISM has a higher density 
in this direction at relatively close distances. 

We find that all the rotational levels can generally be 
fit with a single excitation temperature, T cx ~ 300 K (see 
Table 5), a value that represents all of the sightlines well. 
This follows the results of Spitzer, Cochran, & Hirshfeld 
(1974) where they found if N(J = 0) < 10 15 cm" 2 , a sin- 
gle excitation temperature generally fits all the observed 
N(J). The value of T ex reported in Table 5 compares fa- 
vorably with the results of Spitzer, Cochran, & Hirshfeld 
(1974). Note that they also found one sightline (rSco) 
with a similar total H 2 column density (14.56 dex) but 
with a much larger T ox ~ 1013 K. 

The H 2 b- values for the 6 sightlines in our sample range 
from 2.8 to 5 kms -1 . This is smaller than the 6-values 
derived for O I and N I which are generally larger than 5 
kms _1 (see Table 3). The parameter b is usually defined 
as b — yj2kT/(Am) + w 2 urb , where the symbols have their 
usual meaning. So if H 2 and O I (or N I) were in the same 



7 The distance for WD 1725+586 is from J. Dupuis (2002, private communication). For the other sightlines mentioned in this section the 
distances are from Holberg, Barstow, & Sion (1998) or Vennes et al. (1998) unless otherwise mentioned. 

8 This star was not considered for this survey because the large contamination of molecular lines and stellar lines, but clearly a detailed analysis 
of this sightline would be very interesting. 
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gas phase, one would expect to have 6(0 I) < ^(Eb), since 
the mass of H2 is lower than that of O. This is not the 
case, and it would seem to imply that H 2 is in a differ- 
ent phase. However, some of the higher values for 6(0 I) 
might be explained by the presence of several components 
separated in velocity but not all detected in H2 . The H2 b- 
values indicate a thermal temperature less than 900-3000 
K, different from the temperature of 7000 K observed usu- 
ally in the local interstellar gas. The O I column densities 
and distances of these sightlines place these clouds near 
the wall of the LB or beyond. This would appear to be in 
agreement with H 2 tracing colder neutral gas and hence 
with smaller b- values. 

All the sightlines showing H2 in our survey are located 
in roughly the same area, with similar upper limit dis- 
tances or, as discussed in § 5.2, ranges of distances. They 
also have similar column densities, molecular fractions, ex- 
citation temperatures, and 6-values. These facts suggest 
that the gas studied could be the same large diffuse cloud, 
possibly in form of an extended thin sheet. 

Combining these results imply that the H 2 clouds in the 
LISM can be in form of diffuse clouds, but there are re- 
gions (not explored with our survey) where there is enough 
gas and dust to shield the interstellar radiation field, al- 
lowing a large molecular fraction to form (see, Trapero et 
al. 1996). 

6. IONIZATION OF THE LISM 

6.1. The Neutral Species 
6.1.1. Background 

Recent studies of the neutral species in the LISM reveal 
that they could constrain the origin (photoionization or 
non-equilibrium cooling) of the ionization (Sofia & Jenk- 
ins 1998; Jenkins et al. 2000). The argument of Sofia & 
Jenkins (1998) is that a deficiency of Ar I is not caused by 
Ar depleting into dust grains, but rather its neutral form 
can be more easily photoionized than H. The photonioza- 
tion cross-section of Ar is about 10 times that of H over 
a broad range of energies, so Ar can be more fully ionized 
than H in a partially ionized gas. While in dense regions 
N behaves like O, i.e. the ionization fractions of these ele- 
ments arc strongly coupled with the ionization fraction of 
H via a resonant charge- exchange reaction, in low density, 
partially ionized regions N behaves more like Ar by show- 
ing a deficiency of its neutral form. N also is not depleted 
onto dust grains (Meyer, Cardelli, & Sofia 1997). On the 
one hand, if a region was highly ionized at some earlier 
epoch, by say a shock, and is in the process of recombin- 
ing, then Ar I, (N I) and H I would have about the same 
ionized fraction because the recombination coefficients are 
roughly the same. On the other hand, if steady photoion- 
ization dominates, then Ar I should be deficient from the 
arguments mentioned above. If the region is low density 
and ionized, N I will also be deficient. 

We first describe the results that the data show in this 
section and § 6.2. We will discuss more thoroughly their 
meaning and their implications with new ionization calcu- 
lations in § 6.3. 

Figure 4 presents the ratios of Ar I/O I, N I/O I, and 
(N II+N I)/0 I plotted versus the total column density 
of O I and the distance of the sightline. Both Ar I and 
N I show deficiencies with respect to O I toward many of 



the sightlines but the scatter is large. We also plot in this 
figure the B-type star and solar abundances for Ar/O, and 
B-type star, solar, and interstellar abundances for N/O. 

The reference abundances of Ar are from Kccnan et al. 
(1990) and Holmgren et al. (1990) for B-type stars, and 
from Anders & Grevesse (1989) and Meyer (1989) for the 
solar abundance. The solar O and N abundances are from 
Allende Prieto, Lambert, & Asplund (2001) and Holweger 
(2001), respectively. We note that Holweger (2001) also 
derived an O solar abundance, but with a higher uncer- 
tainty than Allende Prieto, Lambert, & Asplund (2001). 
In any case the difference between these two values is 0.05 
dex and does not change any of our conclusions. The in- 
terstellar O and N abundances are from Meyer, Jura, & 
Cardelli (1998) and Meyer, Cardelli, & Sofia (1997), re- 
spectively (the interstellar abundances were typically mea- 
sured for more distant sightlines with higher hydrogen col- 
umn densities, 20.2 < logA^H I) < 21.2). For the O and 
N abundances in B-type stars see Sofia & Meyer (2001) 
and references therein. Table 6 summarizes the Ar I/O I 
and N I/O I ratios measured in this study, normalized to 
the solar abundances for reference. 

6.1.2. Nitrogen and Oxygen 

The N I/O I ratio versus the distance presented in Fig- 
ure 4 does not show any particular trend, probably primar- 
ily because the interstellar clouds are not resolved and the 
exact distance of the clouds remains unknown. However, 
when considering the sightlines for which the clouds are 
well within the LB (log N(0 I) < 15.60), N I/O I is sys- 
tematically down to —1.00 ±0.07 dex (compared to the in- 
terstellar ratio in denser clouds of —0.66 dex and to the so- 
lar value of -0.76 dex). Once log N(0 I) > 15.60 dex, i.e. 
the gas is near or beyond the edge of the LB, the picture 
is far more complicated, with a scatter of N I/O I between 
—1.1 dex and the interstellar value obtained for sightlines 
with much higher column densities; N I/O I varies over 
0.45 dex in the LISM. 

The bottom panel of Figure 4 shows clearly that mea- 
surements or lower limits of (N II+N I)/0 I are much 
closer to a normal interstellar N/O ratio. The deficiency 
of N is therefore explained by ionization of this element. In 
the LB, the amount of N II can be as large as ^3 times the 
amount of N I. Several limits or values of (N II+N I)/0 I 
are above the interstellar N/O ratio, suggesting that O can 
be partially ionized in significant quantities as well. 

The last column in Table 6 gives the ratio of 
N I/(N I+N II). There is a measured value of this ratio 
toward 5 stars. For the four of them which are well inside 
the LB (logA r (0 I) < 15.60) the measured ratio implies 
that from 40% to 70% of N is ionized, while in the sightline 
located near the edge of the LB (log AT (O I) = 15.78) only 
5% of the N is ionized. The upper limits show, depending 
on the direction, that at least 2% to 75% of N can be ion- 
ized. Combining the upper limits and the measurements 
suggests that the ionization conditions are not uniform in 
the LB (even though N I/O I appears fairly constant). 

6.1.3. Argon and Oxygen 

For most sightlines, the Ar I/O I ratio also shows a de- 
ficiency of Ar I in the LB by ~ —0.3 dex with respect to 
the solar value when logA r (0 I) < 15.60 dex, but with 
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a larger scatter than the scatter observed in N I/O I. 
Ar I/O I is near solar for the WD 1634—573 sightline with 
logiV(0 I) = 15.51 dex. Wood et al. (2002a) using GHRS 
G230M observations suggest that there could be at least 
two clouds along WD 1634—573. The amount of N II is ac- 
tually larger than the amount of N I by a factor 2 toward 
WD 1634-573, even though the N I/O I is the highest 
in the LB sample. It seems therefore plausible that one of 
these clouds is largely ionized, while the other one is mostly 
neutral. This explains why N I/O I and Ar I/O I are near 
solar, since these atoms exist predominantly in the neutral 
cloud. In addition, the super-solar and super-interstellar 
ratios for (N II+N I)/0 I and the low N I/(N II+N I) ratio 
indicate a deficiency of the neutral form of O and N, and 
thus the presence of a substantially ionized cloud along 
this line of sight (this is also confirmed by the strong C III 
line, see § 6.2). Hence, a near solar ratio for Ar I/O I does 
not necessarily imply that there is no ionized cloud along 
the sightline. 

Near the edge of the LB and beyond, the scatter of the 
Ar I/O I becomes larger, ranging from a solar value to 
—0.8 dex subsolar. In the Local Cloud within ~ 5 pc from 
the Sun, Ar II is the dominant ion according to the pho- 
toionization calculation by Slavin & Frisch (2002). The 
scatter in the Ar I/O I ratios is large, either plotted with 
respect to the O I column density or the distance and is 
larger than the scatter of N I/O I. 

The Ar II line can be found at 919.781 A, which is within 
the FUSE spectral range. However, the data in this wave- 
length region are usually of lower quality and this line is 
too weak to be detectable in the spectra studied here. Sofia 
& Jenkins (1998) show that Ar was very unlikely to be de- 
pleted onto dust grains. Therefore we conclude that the 
observed deficiency of Ar I found for almost every sightline 
in our study of the LISM is due to photoionization. 

6.1.4. Nitrogen and Argon 

In Figure 5, we compare the ratios of [Ar I/O I] and 
[N I/O I] normalized to their solar values 9 but only for 
data with la error < ±0.30 dex, i.e. only with the higher 
quality measurements. This plot confirms that N I and 
Ar I are deficient compared to their respective solar val- 
ues. With only a few exceptions, [N I/O I] > [Ar I/O I]. 
This result can be explained by the fact that Ar is more 
easily ionized than N because of its larger ionization cross 
section. There are, however, several data points that sug- 
gest that N I could be as efficiently ionized as Ar I. Figure 3 
of Sofia & Jenkins (1998) suggests that for ionizing pho- 
tons energies E > 35 eV, N could be more easily ionized 
than Ar. But with the constraint from C III/C II mea- 
sured in this study and with new ionization calculations, 
it appears difficult to satisfy this condition in the LISM. 
We will return to a discussion of [Ar I/O I] and [N I/O I] 
in the context of a photoionization model in § 6.3. 

6.2. The Ionized Species 

The singly ionized species measured in this work are in 
the dominant ionization stage in both the neutral and ion- 
ized gas, and therefore cannot tell us directly much about 
the ionization conditions in the LISM. There are, however, 

9 Hereafter, [X/Y] refers to an abundance ratio in logarithmic solar 
[X/Y] = log(X/Y)-log(X/Y) . 



detections of C III along nearly every sightline and one de- 
tection of Fe III toward WD 1800+685. The production 
of C III by photoionization requires photons with energies 
greater than the C II ionization threshold of 24.4 eV. At 
energies above 47.9 eV, C III is ionized to C IV. For the 
production of Fe III, the Fe II threshold is 16.2 eV and the 
Fe III threshold is 30.7 eV. As discussed in § 6.3, C III, in 
particular, is a good tracer of ionized gas. 

Unfortunately, usually both the C II and C III lines are 
saturated, and only a few meaningful limits and measure- 
ments are presented in Table 7. In this Table, we also list 
the values or range of values of C III after correction for 
blending of the interstellar line with the stellar C III line. 
We modeled the synthetic spectra by using either the C III 
excited lines at 1175 A or using C IV available with IUE. 
In some cases, we were only able to estimate a limit and 
hence only a range of possible values. 

Toward WD 1634-573 the value of C III/C II is rela- 
tively high, about —0.63 dex, confirming the presence of a 
largely ionized cloud along this sightline, possibly a H II 
region around the star. For all the other sightlines the 
ratio is less than or of the order of —1 dex; and toward 
WD 1314+293 this ratio is very weak, -2.24 dex. The 
sightlines for which we were able to measure C III/C II lie 
well within the LB. However, Tables 4 and 7 show that 
C III appears to be ubiquitous in the LISM toward our 
sightlines. We note, however, that toward several more 
nearby stars, C III was not detected (Rcdficld ct al. 2002), 
implying that either the ionization conditions change in 
more distant sightlines or that C III might arise partly 
from H II regions near the WDs. While C III/C II ap- 
pears to vary from sightline to sightline, unlike what we 
found for N II/N I, the fraction of C III with respect to 
C II is generally small. The recent ionization calculations 
of Slavin & Frisch (2002) of the Local Cloud within 5 pc 
of the sun show that C III/C II— -1.44 dex (model 17, 
their Table 8), in general agreement with our measure- 
ments. Our observations typically probe denser clouds. 
Possibly a significant fraction of the clouds arc neutral re- 
gions, where C II is also dominant. These observations 
demonstrate that C III, while present, is not a dominant 
ion in the LISM. 

While we generally have a good measurement of 
Fe II, Fe III is either undetected or can be blended 
with photosphcric lines. Only one measurement, 

log[N(Fe m/JV"(Fe II)] 1.2 dex, was obtained toward 

WD 1800+685 which probably lies beyond the wall of the 
LB. This result and the other limits are summarized in 
Table 8. Slavin & Frisch (2002) found with their ioniza- 
tion calculations Fe III/Fe II— —1.40 dex (model 17, their 
Table 8), which, as for C III, shows an agreement with our 
measurements and indicates that Fe II is the dominant ion 
in the LISM. 

6.3. Ionization Calculations and Implications 

To obtain more quantitative insights, we developed the 
photoionization code used by Sofia & Jenkins (1998) and 
Jenkins ct al. (2000). The detailed equations and how they 
are solved are fully described in Sofia & Jenkins (1998). 
But here, in all cases, the formulae included the effects of 

units, X and Y represent the column densities of two elements: 
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charge exchange with neutral hydrogen and the modifica- 
tions of the H and He ionization fractions. We modified 
several parameters in the program that created Figure 2 
in Jenkins et al. (2000) in the following ways: 

(1) We replaced the hot gas interface radiation given by 
Slavin (1989) with the spectrum calculated by Slavin & 
Frisch (1998). This new spectrum shows more details at 
high photon energies, and the physical conditions for the 
Local Cloud seem more realistic. The new spectrum is 
about 5 times more intense than the old one. This slightly 
increases the ionization inside the cloud. We are assuming 
that more distant clouds have properties similar to those 
that we know exist for the Local Cloud. Likewise, we as- 
sume that the radiation field due to the WDs is the same 
everywhere and is that reported for the local neighborhood 
by Vallerga (1998) using EUVE observations. 

(2) The upper limit for the H I absorbing column den- 
sities has been increased to somewhat above 10 19 cm~ 2 , 
on the assumption that some lines of sight could be domi- 
nated by monolithic clouds with this column density. This 
upper limit approximately corresponds to the median col- 
umn density of O I that we found in this survey. In all 
likelihood, the higher column density cases are composed 
of several clouds that have less shielding than a few times 
10 19 cm~ 2 of neutral hydrogen. In addition, the complex 
non-spherical geometry of a cloud may reduce the effective 
shielding to absorbing column densities much lower than 
the value implied by the measurement along the sightlinc. 

(3) We added the atomic parameters for C, P, Si, S, and 
Fe. In particular, we plot explicitly the values of C III/C II 
and Fe III/Fe II. 

The results are plotted in Figure 6 as a function of the 
absorbing column shielding in the external ionizing radi- 
ation sources. A uniform pressure within a given cloud 
of p/k — 2.5 x 10 3 cm~ 3 K is assumed. This pressure 
is consistent with the pressure derived by Gry & Jenkins 
(2001) toward eCMa. Reducing the pressure would give 
a C III/C II ratio that is too high with respect to the 
observed values. 

Figure 6 illustrates how O I follows H I very well, while 
N I and Ar I become progressively more deficient as the 
shielding from hydrogen decreases. The other elements 
(C, P, Si, S, and Fe) remain primarily in the singly ionized 
state, but the abundances of these ions are elevated rel- 
ative to that of neutral hydrogen. The assumed pressure 
seems to give a reasonable set of values for C III/C II, but 
Fe III/Fe II appears to be too low. Yet, we have only one 
measurement and this measurement is very uncertain (see 
§ 6.2). In contrast, Slavin & Frisch (2002) found similar 
fractions for C III/C II and Fe III/Fe II in their model. 
Except for this discrepancy, these two models are in rea- 
sonable agreement. According to our model, the fraction 
of Fe III is heavily dependent on the effect of charge ex- 
change with neutral hydrogen. This is not the case for 
C III. Without taking into account the effect of charge ex- 
change with neutral hydrogen, the fraction of Fe III is com- 
parable to the fraction of C III. Therefore, the amount of 
Fe III depends on the density and ionizing radiation field, 
making Fe III a less reliable tracer of ionized gas compared 
to C III. 

In Figure 5, the model predictions are compared with 
the measured ratios of [N I/O I] and [Ar I/O I]. There is 



a large scatter in the measurements due to the scatter in 
the physical properties of the different clouds. The model 
assumes the same physical properties for all clouds. The 
study reported here and many others (e.g., Gry & Jenkins 
2001; Rcdficld et al. 2002) demonstrate that this is not 
the case. Nevertheless, it should be noted that, within 
la error, the majority of the data points agree with our 
model. 

For about 6/22 data points in Figure 5, [Ar I/O I] is 
significantly higher than predicted and thus is approxima- 
tively equal to [N I/O I]. This suggests that, for these 
sightlines, the ionization fraction of N could be as large 
as that of Ar. Using our model, we have explored sev- 
eral combinations of parameters by varying the pressure 
and ionizing flux to see if we could reproduce such ion- 
ization fractions. First, the value of p/k must be well be- 
low 10 3 cm -3 K to make the electron density high enough 
(and n(H I) low enough) in order to have the charge ex- 
change relatively inefficient between nitrogen and hydro- 
gen. Moreover, the spectrum of the ionizing radiation 
must be altered drastically, e.g., with a very strong flux of 
photons (of unknown origin) that have an energy E > 35 
cV (see Figure 3 of Sofia & Jenkins 1998). Perhaps these 
photons originate from the recombination of He II. When 
all of these extreme conditions are invoked, the ionization 
of N can even be slightly higher than that of Ar. How- 
ever, under those conditions C III/C II becomes far too 
large compared to our result. So with our model, a higher 
or comparable ionization fraction of N compared to that 
of Ar is difficult to explain with a simple photoionization 
phenomena. Assuming collisional ionization equilibrium 
in a hot plasma, Shull & van Steenberg (1982) calculated 
the expected deficiency of the neutral forms of Ar, N, and 
O. According to their values, collisional ionization of these 
elements would not appear to solve this problem. While 
in principle O can be depleted by about 25% (Moos et al. 
2002), Ar and N are very unlikely to be depicted by any 
significant amount (Sofia & Jenkins 1998; Meyer, Cardclli, 
& Sofia 1997). Another possibility would be an overabun- 
dance of Ar in certain directions. But that would seem to 
be difficult to reconcile with no apparent changes in abun- 
dances for the other elements, especially the a-elements 
such as O. 

While our model might be too naive to reproduce the 
different changes in the physical conditions occurring in 
the LISM, our observations strongly indicate that pho- 
toionization is a major ionization contribution in the 
LISM. These results do not support the proposition that 
the LISM is in the process of recombining from an highly 
ionized phase at some earlier epoch, in which case Ar I 
and N I and H I would be roughly equally ionized. The 
photoionization conditions may not be constant in the LB 
and more generally in the LISM according to the vari- 
ation observed in Ar I/O I, N I/O I, N I/(N I+N II), 
and C III/C II. But this is complicated by the unresolved 
cloud structure. For example, toward WD 1634—573, at 
least two clouds are present with one of them being sub- 
stantially ionized (large C III and N II fractions; possibly 
an H II region near the star) and the other one mostly 
neutral (near solar abundance for Ar I/O I and N I/O I). 
Another example is toward WD 1314+293, where there is 
a large fraction of N II but a very small fraction of C III. 
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7. DUST IN THE LISM 

Several studies have examined the deficiency of elements 
in the gas-phase of the Galaxy with respect to solar or 
B-type star abundances (see the review by Savage & Sem- 
bach 1996a, and references therein). When ionization is 
not the cause, these deficiencies are called depletions and 
are usually explained by the fact that certain elements are 
preferentially locked up into dust grains. These studies 
show a general progression of increasingly severe depletion 
from warm halo clouds, to warm disk clouds, to colder disk 
clouds. Below, we will compare these observed depletion 
patterns with studies of other parts of the Galaxy. S II 
(or Zn II) or H I are often used because the former two 
are lightly or not depleted onto dust grains and the lat- 
ter gives the absolute gas-phase abundance directly (when 
ionization is not important). S II and Zn II are not avail- 
able in the FUSE wavelength range, and as discussed in 
the § 6, the LB and LISM are partially ionized, there- 
fore using O I as a proxy for H is not suitable. However, 
we show in § 6.3 that as long as conditions are not ex- 
treme, the other elements (C, Si, P, and Fe) appear to 
remain primarily in the singly ionized state. This is also 
confirmed by the photoionization calculations of Slavin & 
Frisch (2002). Therefore the abundance of C II, Si II, P II, 
and Fe II can be compared directly in order to measure 
the depletion onto grains. We note that this method of 
sensing the existence of dust relies on the underlying as- 
sumption that the reference abundances of the elements 
apply to the LISM. Direct detection of dust grains within 
the solar system were made possible using the Ulysses and 
Galileo satellites and are discussed by Frisch et al. (1999). 

Previous studies show that P is only lightly depleted 
onto dust grains especially in warm gas (e.g., Welty et 
al. 1999). However, P II is 2 or 3 orders of magnitude 
less abundant than the other species, so that only a few 
measurements were possible in our survey. On the other 
hand, C is very abundant and is lightly depleted. Because 
of this and the fact that C is so abundant and the only 
available transition of C II is so strong, a measurement of 
the column density of C II was possible toward only a few 
sightlines with low columns (the high velocity components 
have all low column densities and are discussed in § 9). 
Si II and Fe II are both depleted into dust grains, and 
they are easily measured in our survey. 

Evidence for dust grains in the LB and beyond is 
shown in Figure 7 where the logarithmic ratios of 
N(Fe II)/7V(Si II) and N(P II)/iV(Si II) are shown. Even 
though P is the least depleted amongst these three ele- 
ments and would then be a better element for reference, 
we chose Si because there are many more measurements 
in common with Fe. The solar (meteoric) values of these 
ratios are —0.04 and —1.98 dex for Fe/Si and P/Si, respec- 
tively (Anders & Grevesse 1989). Except for one value that 
could be compatible with a solar value within the error 
bars (WD 2011+395), [Fe Il/Si II] is depleted from about 
—0.3 to —0.8 dex, typical of what is observed in "cold- 
warm" and "halo" clouds (Savage & Sembach 1996a, and 
references therein). Note that for Fe/Si it is not possible 
to differentiate between a cold and a warm depletion pat- 
tern, since it is depleted in both types of gas by ~ — 1 dex. 
Figure 7 shows that there is no difference between the in- 
terstellar clouds within the LB, near the wall, or beyond. 



In particular the range of values over which [Fe II /Si II] 
varies is similar. 

[P Il/Si II] varies generally from a solar value (—2.0 dex) 
to —0.8 dex. This is also similar to depletions observed ei- 
ther in "cold" or "warm-halo" clouds (Savage & Sembach 
1996a). Note that in this case the warm and halo depletion 
patterns can not be differentiated from each other (there 
is only an estimate of the depletion of P in halo gas). Only 
two values are available within the LB: a solar value and 
—0.3 dex below solar. 

Plotting the results for [P Il/Si II] and [Fe Il/Si II] for 
common sightlines could a priori differentiate between the 
different depletion patterns since the cold, warm, and halo 
depletion patterns have different locii. Figure 8 shows 
[P Il/Si II] versus [Fe Il/Si II] along with the expected 
depletion patterns observed in more distant sightlines in 
the Galaxy. Most of the clouds appear to be dominated by 
more warm-like or halo-like gas than cold gas. The scatter 
is probably mainly explained by the fact that the clouds 
are not resolved and therefore are a combination of dif- 
ferent gas phases. The sightline toward WD 0131 — 163 is 
more complicated to understand because [P II /Si II] points 
toward a cold depletion pattern while [Fe Il/Si II] toward 
an halo depletion pattern. 

The solar ratios of Si/C, P/C, Fe/C are 
-1.00, -2.98, -1.04, respectively. Toward the few sight- 
lines where a measurement of C II was possible, the ra- 
tios ([Si II/C II], [P II/C II], [Fe II/C II]) are: (0.18, < 
1.40, < -0.13) toward WD 0549+158, (< 0.05, < 1.64, < 
-0.06) toward WD 1254+224, (-0.98, < 0.12,-1.62) 
toward WD 1314+293, (-1.14,-0.84,-1.78) toward 
WD 1634—573. The two first sightlines do not have any 
stringent limits. Toward WD 1314+293, these ratios are 
compatible with a cloud having a cold depletion pattern 
(in agreement with the ratios of P/Si and Fe/Si). These 
ratios for WD 1634—573 point also to a cold depletion 
pattern ([P II/C II] is, however, far too small), but this is 
in contradiction with the results from P/Si and Fe/Si (see 
Figure 8), which might suggest that there is a lack of dust 
containing C toward this sightline. 

The main reason for the deficiencies of Si and Fe are 
understood to be due to the incorporation of some of the 
Fe and Si into dust grains. The scatter in the values may 
reflect changes in the dust fraction or complications in the 
ionization fractions that are similar to those that we found 
for N, Ar, and O. The fact that we do not observe any trend 
with the distance or total column density of the sightline 
is probably mainly due to the fact that the clouds are not 
resolved. 

8. ELECTRON DENSITY AND COOLING RATE IN THE LISM 

8.1. Electron Density 

If collisions with electrons are the principal source of 
fine-structure excitation of C II in the LISM, then the fa- 
miliar collisional equilibrium equation between C II and 
C II* (e.g., Gry & Jenkins 2001) can be simplified to 

n e = 0.189 T°-«^ cm- (1) 

when N(C+*) <C N(C+) and T > 93 K. This relation pro- 
vides an estimate of the electron density when the space 
densities are replaced by column densities. The collisions 
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with neutral hydrogen atoms should not be significant in 
the LISM. Indeed, according to Fig. 3 of Keenan et al. 
(1986), n(H I) = 1.0 cm -3 alters the inferred electron 
density only when log[A(C+*)/A(C+)] < -2.7, but our 
calculations indicate that n(H I) < 0.3 cm -3 in the LISM 
as long as p/k < 2500 cm -3 K (see Figure 6). 

There is no evidence of C being depleted into dust in 
the LISM (see § 7). P is also known to be only lightly de- 
pleted. Moreover, in other regions of the Galaxy, C and P 
have similar depletions (Savage & Sembach 1996a). Differ- 
ences in ionization of C and P are probably not important 
(N(C III) is small with respect to N(C II); see § 6 and 
Figure 6), so that we can write log[iV(C II*)/iV(C II)] « 
logiV(C II*) - log7V(P II) + log(P/C) , where the solar 
ratio log(P/C) Q = -3.0 dex. 

Under these assumptions and assuming a temperature 
of T = 7000 K observed in LIC, we can determine the 
electron density. This is shown in the bottom panel of 
Figure 9. The mean value is n e < 0.1 cm -3 . Toward 
Capella situated at a distance of d = 12.9 pc, Wood & 
Linsky (1997) determined n e = 0.05-0.23 cm" 3 . We also 
plot in Figure 9 the expected electron density from our 
model versus the shielding of H I (see § 6) . All the points 
fall above this line, suggesting that they are probably mul- 
tiple regions where the shielding in each region is much less 
than the total column density of the sightlinc. 

8.2. Cooling Rate 

Because carbon has a high abundance, one of the most 
efficient way for cooling the diffuse gas is via the excitation 
of the C II by collisions with electrons and hydrogen atoms 
followed by the spontaneous emission at 157.7 fim when 
the levels decay to the ground state. It follows that the 
C II cooling rate in ergs -1 per H I atom can be written 
(e.g., Savage & Sembach 1996b), 

lo = 2.89 x 10-"^ , (2) 

where log7V(H i) = \ og N{0 I) + 3.5 (sec § 4). In 
those conditions, the mean value in our sample is log/ c = 
-25.59 ± 0.32 dex. Using P II as a proxy for the to- 
tal H (H I+H II) and assuming no depletion, we have 
logic = —25.30 ± 0.39 dex in our sample. If a depletion 
of ~ —0.2 to —0.3 dex is assumed for P (Savage & Sem- 
bach 1996a; Welty et al. 1999), both determinations are 
in good agreement. This result is in agreement with pre- 
vious determinations derived in the diffuse gas in other 
parts of the Galaxy via different methods (Gry, Lequeux, 
& Boulangcr 1992; Bock et al. 1993; Savage & Sembach 
1996b; Boulanger et al. 1996). Toward Capella, using the 
H I and C II* column densities from Wood et al. (2002b), 
logl c = —25. 14 ±0.03 dex, but increases by about 0.3 dex 
when using O I as a proxy of H I. 

In the top panel of Figure 9, we show the cooling rate 
against the total H I column density, showing a clear anti- 
correlation between these two quantities. This likely indi- 
cates that for low column density sightlines, a substantial 
amount of C II is associated with ionized gas, which is 
not taken into account in the normalization by the N(H) 
column density. 

9. HIGH VELOCITY COMPONENTS 



The sightlines toward WD 0455-282, WD 1202+608, 
and WD 1528+487 each show a high velocity component 
separated by —56, —40, and —39 kms -1 , respectively, with 
respect to the lower velocity component. In the different 
tables, the high velocity component is designated by the 
letter "a" after the star name. Toward WD 0455— 282 and 
WD 1528+487, high velocity components are clearly de- 
tected only in the C II absorption line for the elements 
considered here. For WD 0455— 282, there is also a clear 
detection of this component in the H I Lyman lines. For 
WD 1528+487, there is only a hint in the Lyman lines. 
This rules out misidentification due to a stellar line. For 
these two stars, the presence of C II absorption features 
and the non-detection of neutral species, in particular the 
strong O I A1039 line, implies that these clouds are mostly 
ionized. 

Toward WD 1202+608, several species are detected (see 
Table 4) , making the identification as a high- velocity com- 
ponent unquestionable. In contrast to the two other sight- 
lines which have a very low column density with respect to 
their low velocity counterparts, this cloud has comparable 
column densities compared to the lower velocity compo- 
nents for Si II and Fe II. However, it has much smaller 
column densities for the neutral species, O I and N I, im- 
plying that this cloud is predominantly ionized. In addi- 
tion, N I/(N II+N I) < 0.09, showing that at least 91% 
of N is in an ionized form. A similar comment applies 
for Ar I, even though the result remains more uncertain 
for this atom. Unfortunately, we have no information for 
C III and Fe III, because the former is blended with the 
local component (but the large lower limit - the largest in 
our sample - suggests a fair amount of C III) , and Fe III is 
confused by the moving stellar lines from this binary sys- 
tem. The column density of O I suggests an H I column 
density of ~ 18.20 dex, while Si II gives a total H column 
density of > 18.6 dex (it is a lower limit because Si can be 
substantially depleted onto dust grains), implying an H II 
column density > 18.3 dex. 

Holberg et al. (1995) and Kruk, Chayer, & Dupuis 
(2003) argue that the origin of the high velocity component 
toward WD 1202+608 may be a remnant from the ejected 
common envelope of the binary, assuming it is a relatively 
young system. Gry & Jenkins (2001) also identified two 
negative high velocity components in the line of sight to- 
ward eCMa at —27 kms -1 and —82 kms -1 from the main 
component. The fact that the velocity of these components 
is always negative would also point toward a circumstellar 
origin. However, it should be noted that the high velocity 
features appear in pairs of stars. WD 0455— 282 is in the 
same part of the sky as e CMa, and WD 1202+608 is in 
the same part of the sky as WD 1528+487. This suggests 
that in each of these two cases, the sightlines could be 
passing through a single gas cloud. If so, the clouds are 
probably much closer than the stars and because of the 
high ionization are likely well within the LB. Although the 
argument is not definitive, it does weigh against the hy- 
pothesis of stellar ejecta. Finally, we note that this study 
does not show widespread high- velocity components along 
many sight lines. In particular, there is no evidence of 
globally expanding gas at positive velocities as might be 
expected for the case of a recent supernova explosion po- 
sitioned near the Sun. 
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10. SUMMARY AND CONCLUDING REMARKS 

We present a comprehensive survey of the LISM which 
includes the main species available in the far-uv wave- 
length range (C II, C II*, C III, N I, N II, O I, Ar I, 
Si II, P II, Fe II, Fe III, and H 2 ). The spectral resolution 
(R ~ 20, 000) is not sufficient for detailed analysis of veloc- 
ity structures. However, the number of species available 
permit a study of some of the fundamental characteris- 
tics of the local interstellar gas by comparing the relative 
column densities, the ionization structure, the molecular 
content, and the indirect indication of dust. 

Our survey detected only very diffuse molecular clouds 
(/h 2 < 10~ 4 ). H2 was detected in only six out of 31 sight- 
lines all within a single sector of the sky (see Figure 1). 
Thus it appears that H 2 is not widely distributed in the 
LISM. The derived H2 6-values for the 6 sightlines in our 
sample range from 2.8 to 5 kms -1 , implying thermal tem- 
peratures less than 900-3000 K. This is different from the 
temperature of 7000 K usually observed in local interstellar 
clouds, implying different gas phases in the LISM. Com- 
bining our results with previous work on H2 in the LISM 
shows that the molecular phases are not uniform. Molecu- 
lar hydrogen can exist in diffuse clouds, but in some regions 
the molecular fractions are similar to the atomic content. 

The comparison of column densities of Si II, P II, Fe II 
provides information on the dust depletions of LISM. The 
relative abundances of these species vary from sightlinc 
to sightline, again suggesting different gas-phases in the 
LISM. These abundance ratios are better represented by 
typical warm or halo diffuse clouds observed in more dis- 
tant sightlines in the Galaxy rather than a cold depletion 
pattern. 

The cooling rate in the LISM, l c (in ergs -1 H I atom -1 , 
where the H I column density was determined from the 
O I column density), can be derived directly from the 
C II* absorption lines and has a mean value of log^ c = 
— 25.59±0.32 dex, very similar to previous determinations. 
Higher values are observed at lower H I column densities, 
indicating that for low column density sightlines, a sub- 
stantial amount of the C II must come from ionized gas. 
The electron density can also be inferred from the C II* 
absorption lines. Using P II as a proxy of C II and assum- 
ing a temperature of 7000 K, the electron density in the 
LISM is n e < 0.1 cm -3 , but larger than the values pre- 
dicted by our model. This disagreement likely indicates 
that there are multiple regions along the sightlines where 
the H I shielding in each region is much less than the total 
H I column density of the sightline. 

The comparison of column densities of N I and Ar I 
to O I provides information on the ionization structure 
of the LISM. The N I/O I ratio is systematically -0.24 
dex subsolar in the LB. Near or beyond the LB bound- 
ary, the picture is far more complicated, with a scatter 
of N I/O I between ~ —0.3 dex subsolar to ~ +0.15 dex 
supersolar (WD 1636+351). Using measured values and 
lower limits on the N II column densities, we show that 
the observed deficiencies are likely due to ionization and 
not abundance variations or dust depletions. For most 
sightlines, the Ar I/O I ratios also show deficiencies of 
Ar I in the LB by ~ —0.35 dex with respect to the solar 
value, but with a larger scatter than for N I, and this scat- 
ter increases beyond the LB ranging from a solar value to 



~ — 0.8 dex subsolar. 

More highly ionized species such as C III and Fe III are 
also available to constrain the ionization structure. C III 
appears to be ubiquitous in the LISM for our sample, and 
C III/C II appears to vary from sightline to sightline. The 
C III/C II ratio is however generally small; C III is not 
a dominant ion in the LISM. Fe III remains more elusive 
and is often completely lost in the stellar profile. But our 
limits on Fe III/Fe II indicate that Fe II is the dominant 
ion in the LISM. 

We modified the photoionization code used by Sofia & 
Jenkins (1998) and Jenkins et al. (2000). In particular, 
we adopted the new radiation field calculated by Slavin & 
Frisch (1998) and added the atomic parameters for C, P, 
Si, S, and Fe. The calculations were performed for a se- 
quence of models where the thermal pressure was set to a 
value of p/k = 2.5 x 10 3 cm -3 K. These calculations show 
that O I follows H I very well, while N I and Ar I become 
progressively more deficient as the shielding from H I be- 
comes smaller. The other elements (C, P, Si, S, and Fe) 
remain primarily in the singly ionized state. The assumed 
pressure also seems to give a reasonable set of values for 
C III/C II. However, the values calculated for Fe III/Fe II 
appear to be too low, although we have only one measured 
value of Fe III/Fe II, the rest are upper limits. According 
to our model, the fraction of Fe III is heavily dependent on 
the effect of charge exchange with neutral hydrogen and is 
therefore not a reliable tracer of ionized gas. 

Ar I, N I and H I have similar recombination coefficients, 
so if the LISM was highly ionized at an earlier epoch by 
a strong shock and is in the process of recombining, these 
elements would have roughly equal ionized fractions. Our 
study contradicts this proposition by indicating that pho- 
toionization is a major contribution. We note, however, 
that in a few cases, [N I/O I] is approximatively equal 
to [Ar I/O I], suggesting that the ionization fraction of 
N could be as large as that of Ar, which contradicts our 
model. This could indicate the radiation field might be 
harder in some directions. However, our model does not 
provide a clear explanation for this phenomenon with re- 
spect to the different observational constraints, and other 
factors are needed to explain this condition (e.g., non- 
equilibrium collisional ionization or overabundance of Ar 
possibly due to a supernova). 

The variations observed in the ratios of Ar I/O I, 
N I/O I, N I/(N I+N II), and C III/C II indicate that 
the photoionization conditions are probably not uniform 
in the LB and in the LISM. Such changes in the pho- 
toionization conditions might be linked to the patchincss 
observed in the distribution O VI absorption (Oegerle et 
al. 2003). These variations, in turn, may depend on the 
strength and direction of the magnetic field at the cloud 
boundaries. Slavin & Frisch (2002) argued that evapora- 
tive interfaces at the boundaries of the clouds (O VI ab- 
sorption being a tracer of such interfaces) could account 
for a large fraction of the observed ionization. Better ob- 
servational constraints will come from higher spectral res- 
olution of far-UV spectra, in which the individual clouds 
can be resolved. 

Finally, we note that our results do not strongly con- 
tradict the X-ray data. Juda et al. (1991) compared the 
Be and B band fluxes in a number of different directions 
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and found that the ratio of the two did not change with 
the overall intensity. From this nearly linear relationship 
between Be and B band fluxes, they concluded that inter- 
spersed neutral gas within the LB must have an average 
density less than 0.025 cm~ 3 . Unfortunately, their sky cov- 
erage was limited and most of their viewpoints were not 
near any of our targets. An exception is WD 0501+527, 
for which logiV(H I) = 18.11 (Lemoine et al. 2002), and 
d = 69 pc, which gives (n(R I)) = 0.007 cm~ 3 , below 
their limits. But only 30° away from one of the X-ray di- 
rections, toward WD 1631+781, the same calculation gives 
0.1 cm~ 3 . It would be interesting to see if toward denser 
lines of sight (e.g., WD 0004+330 gives (n(H I)) = 0.24 
cm -3 ), there is measurable X-ray absorption in the lowest 
energy band. 
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Table 1 

SlGHTLINE SUMMARY 



WD name 


Other name 


; 

(°) 


b 

(°) 


fholio 

(kms -1 ) 


d 

(PC) 


Spectral 
Type 


Aperture 


Notes 


WD 0004+330 


GD 2 


112.48 


-28.69 


+0.1 


97 a 


DA 


LWRS, MDRS 


A 


(1) 


WD 0050-332 


GD 659 


299.15 


-84.12 


+ 12.2 


58 a 


DAI 


LWRS 


B 


(2) 


WD0131-163 


GD 984 


167.26 


-75.15 




96 a 


DAI 


LWRS 


B 


(2) 


WD 0455-282 


MOT 0455-2812 


229.29 


-36.17 


+ 13.6 


102 a 


DAI 


MDRS 


C 


(2) 


WD 0501+527 


G191-B2B 


155.95 


+7.10 


+ 13.0 


69 ± 16 (59°) 


DAI 


LWRS, MDRS, HIRS 


C 


(3) 


WD 0549+158 


GD 71 


192.03 


-5.34 


+23.2 


49 a 


DAI 


LWRS 


A 


(2) 


WD 0621-376 




245.41 


-21.43 


+ 19.5 


78 a 


DA 


LWRS, MDRS 


D 


(4) 


WD 0715-703 




281.62 


-23.50 




94" 


DAI 


MDRS 


A 


(2) 


WD 1017-138 




255.74 


+34.53 


+34.5 


90 a 


DA 


LWRS 


A 


(2) 


WD 1202+608 


GD 314, Feige 55 


133.11 


+55.66 


-16.1 


200 a 


DA 


MDRS 


C 


(2) 


WD 1211+332 


HZ 21 


175.03 


+80.02 


-18.0 


115 ± 35 


DO 


LWRS, MDRS 


B 


(5) 


WD 1234+481 




129.81 


+69.01 


-28.9 


129° 


DA 


LWRS 


A 


(2) 


WD 1254+223 


GD 153 


317.26 


+84.75 


-11.1 


67" 


DAI 


LWRS 


A 


(2) 


WD 1314+293 


HZ 43 


54.10 


+84.16 


-17.4 


68 ± 13 


DAI 


LWRS, MDRS 


A 


(6) 


WD 1528+487 




78.87 


+52.72 




140° 


DA: 


LWRS 


B 


(2) 


WD 1615-154 


EGGR 118 


358.79 


+24.18 


-38.2 


55 a 


DAI. 5 


MDRS 


A 


(2) 


WD 1631+781 


1ES 1631+78.1 


111.29 


+33.58 


-11.8 


67 a 


DAI 


MDRS 


A 


(7) 


WD 1634-573 


HD 149499 B 


329.88 


-7.02 


-19.6° 


37 ±3 


DOZ1 


LWRS, MDRS 


B 


(8) 


WD 1636+351 




56.98 


+41.40 




109° 


DA 


LWRS 


A 


(2) 


WD 1800+685 




98.73 


+29.78 


-15.9 


159° 


DAI 


MDRS 


A 


(2) 


WD 1844-223 




12.50 


-9.25 




62 a 


DA: 


LWRS 


A 


(2) 


WD 2004-605 




336.58 


-32.86 


-28.0 


58° 


DA 


MDRS 


A 


(7) 


WD 2011+395 


EUVE J2013+40.0 


77.00 


+3.18 




141 a 


DAO 


LWRS 


B 


(2) 


WD 2111+498 


GD 394 


91.37 


+ 1.13 


-6.2 


50 a 


DA2 


LWRS 


B 


(2) 


WD 2124-224 




26.81 


-43.19 




224 6 


DA 


MDRS 


B 


(2) 


WD 2148+286* 


BD+28°4211 


81.87 


-19.29 


-20.2' J 


104 ± 18 


SdO 


LWRS, MDRS 


E 


(9) 


WD 2152-548 




339.73 


-48.06 




128° 


DA: 


MDRS 


A 


(2) 


WD 2211-495 




345.79 


-52.62 


-0.7 


53" 


DA 


LWRS, MDRS, HIRS [Dl (10) 


WD 2247+583 


Lan 23 


107.64 


-0.64 




122° 


DA 


1WRS 


C 


(5) 


WD 2309+105 


GD 246 


87.26 


-45.11 


-7.9 


79" 


DAI 


LWRS, MDRS 


E 


(5) 


WD 2331-475 


MCT 2331-4731 


334.83 


-64.81 


+ 14.9 


82 a 


DA 


1WRS, MDRS 


B 


(5) 



Note. — Heliocentric velocities arc from Holbcrg, Barstow, & Sion (1998), except (a) from Wood et al. (2002a) and (/3) from Sonneborn 
et al. (2002). *Notc that using the denomination WD 2148+286 for BD-f28°4211 is somewhat incorrect because this star is a subdwarf 
of type O, not a white dwarf. However, for consistency and simplicity, we have used this name throughout the article, a: Distances from 
Vennes et al. (1997). b: Distance from Vennes et al. (1998). c: Distance from Dupuis (2002, private communication), d: Distances from 
Holbcrg, Barstow, & Sion (1998). Distances with errors arc from Hipparcos parallax measurements (Pcrryman et al. 1997); for the others, 
the errors are typically 20—30%. Metallicity of the star: [A] no metal detected in FUSE spectra; [B] few metals; [C] rich in metals; [D] 
very rich in metals; [E] solar metallicity (sec Dupuis et al. 2003, in preparation, and references below). Source of the column densities: (1) 
Olivcira et al. (2003, private communication), this work; (2) this work; (3) Lcmoine et al. (2002); Andre et al. (2003b); (4) Lehner et al. 
(2002), this work; (5) Oliveira et al. (2003), this work; (6) Kruk et al. (2002); (7) Hebrard et al. (2003), this work; (8) Wood et al. (2002a), 
this work; (9) Sonneborn et al. (2002), this work; (10) Hebrard et al. (2003), this work. 
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Table 2 

Atomic data used in this work 



Ions 


(A) 




/ 






Note 


ClI 


1036.337 


1.18 


X 


10" 


-1 


saturated (a) 


ClI* 


1037.018 


1.18 


X 


10" 


-1 




c III 


977.020 


7.58 


X 


10" 


-1 


blended (b) 


Nl 


952.303 


1.70 


X 


10" 


-3 






952.415 


1.87 


X 


10" 


-3 






953.415 


1.32 


X 


10" 


-2 






953.655 


2.50 


X 


10" 


-2 






953.970 


3.48 


X 


10" 


-2 






963.990 


1.48 


X 


10" 


-2 






964.626 


9.43 


X 


10" 


-3 






965.041 


4.02 


X 


10" 


-3 






1134.165 


1.52 


X 


10" 


-2 


airglow (c) 




1134.415 


2.97 


X 


10" 


-2 


airglow (c) 




1134.980 


4.35 


X 


10" 


-2 


airglow (c) 


Nil 


1083.994 


1.15 


X 


10" 


-1 


saturated (a) 


Ol 


919.658 


9.48 


X 


10" 


-4 






919.917 


1.78 


X 


10" 


-4 






921.857 


1.19 


X 


10" 


-3 








2.45 


X 


10 


-4 






924.950 


1.54 


X 


10" 


-3 






925.446 


3.50 


X 


10" 


-4 






929.517 


2.29 


X 


10" 


-3 






930.257 


5.37 


X 


10" 


-4 






936.630 


3.65 


X 


10" 


-3 






948.686 


6.31 


X 


10" 


-3 






950.885 


1.58 


X 


10" 


-3 






976.448 


3.31 


X 


10" 


-3 






1039.230 


9.20 


X 


10" 


-3 


airglow (c) 


Si II 


1020.699 


1.64 


X 


10" 


-2 




Ar I 


1048.220 


2.63 


X 


10" 


-1 






1066.660 


6.75 


X 


10" 


-2 


blended (d) 


PlI 


1152.818 


2.45 


X 


10" 


-1 




Fell 


1055.262 


7.50 


X 


10" 


-3 






1063.176 


5.47 


X 


10" 


-2 






1096.877 


3.20 


X 


10" 


-2 






1125.448 


1.60 


X 


10" 


-2 






1143.226 


1.77 


X 


10" 


-2 






1144.938 


1.06 


X 


10" 


-1 




FoIII 


1122.524 


5.44 


X 


10" 


-2 





Note. — Rest frame vacuum wavelengths and os- 
cillator strengths are from Morton (private communi- 
cation, 2000), except for the /-values of Fe II which 
are from Howk et al. (2000). (a) These lines are gen- 
erally strong and usually only a lower limit can be 
estimated (except in few cases where the column den- 
sities are extremely small), (b) Interstellar C III can 
be blended with the stellar component if the star con- 
tains some metals, (c) Those lines can be affected by 
airglow emission lines, especially the N I triplet, (d) 
Ar I A1066 can be blended with stellar Si IV A1066 if 
the star contains some metals. 
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Table 3 

Comparison of O I and N I 6-values 



Sightline 



b [OI] 
(kms 1 ) 



b [N I] 
(kms^ 1 ) 



WD 0004+330 
WD 0050-332 
WD 0131-163 
WD 0455-282 
WD 0715-703 
WD 1017-138 
WD 1202+608 
WD 1234+481 
WD 1528+487 
WD 1615-154 
WD 1631+781 
WD 1636+351 
WD 1800+685 
WD 1844-223 
WD 2004-605 
WD 2111+498 
WD 2124-224 
WD 2152-548 
WD 2331-475 



6.3+ ° 

5.3+ » 

5.5+ » 

5.8+ « 

6.4+ » 

5.7+ a 

5.4 + 

5.o± a 
5.8+ a 

4.3 : 
6.4+ ° 
6.1 + 

8.4 + 

4.5 + 
4.3 + 
8.0+ 10 
10.7+ oj 
5.0 + 

8.6 + 



hi 

8:1 

8:3 

n 



0.4 
0.3 



M 

0.4 



4.8 + 
5.3 
5.5 

6.2 ± 

6.2 ± 
4.8 + 

5.0 

6.3 + 
4.3 + 
4.7 + 

5.7 + 

8.8 + 
4.3 + 
4.5 ± 

7.9 + 
10.7+ i 
5.4+ 0; 
7.8+ k 



Note. — ":" means that the fe-value was 
taken from the other species; see text. 
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Table 4 
Column density summary 



Star 


ClI 


ClI* 


ClII 


Ol 


Nl 


Nil 


Sill 


Ar I 


Foil 


FoIII 


PlI 


WD 0004+330 


>14.38 


>13.51 


>13.40 


16.20 
16.35 
16.50 


15.15 
15.30 
15.45 


>14.12 


> 14.34 


>13.80 


13.55 
13.60 
13.68 


<12.95 


12.83 
12.85 
12.87 


WD 0050-332 


>14.09 


12.29 
12.55 
12.71 


<13.14 


15.21 
15.26 
15.32 


14.01 
14.30 
14.60 


>13.91 


<13.60 


<12.67 


12.80 
13.06 
13.22 




<12.39 


WD 0131-163 


>14.20 


12.84 
13.00 
13.12 


<13.27 


15.83 
15.89 
15.97 


14.80 
14.91 
15.04 


>14.11 


13.80 
13.93 
14.05 


12.90 
12.97 
13.03 


13.51 
13.54 
13.57 




12.66 
12.78 
12.87 


WD 0455-282a 


13.54 
13.58 
13.62 


<12.68 


<12.66 


<14.12 


<13.13 


<13.47 


<12.88 


<12.62 


<12.73 




<12.39 


WD 0455-282b 


>14.20 


12.46 
12.68 
12.83 


>13.72 


14.91: 


13.29 
13.44 
13.59 


>13.92 


13.34: 


<12.62 


12.83 
12.94 
13.02 


<12.99 


<12.39 


WD 0501+527 


>15.39 


13.23 
13.26 
13.28 




14.82 
14.86 
14.90 


13.83 
13.87 
13.91 


>13.83 


13.51 
13.52 
13.53 


<12.00 


13.03 
13.05 
13.07 




<11.85 


WD 0549+158 


13.79 
13.81 
13.84 


<12.84 


12.59 
12.72 
12.81 


14.18 
14.27 
14.34 


<13.60 


13.43 
13.49 
13.54 


12.75 

12.99 
13.14 


<12.46 


<12.64 


<13.12 


<12.23 


WD 0621-376 


>14.19 


13.04 
13.09 
13.14 


<13.31 


15.22 
15.26 
15.31 


14.27 
14.34 
14.41 


>13.80 


<13.68 


12.79 
12.86 
12.93 


12.79 
12.94 
13.05 






WD 0715-703 


>14.29 


13.31 
13.42 
13.51 


>13.22 


15.81 
15.90 
16.00 


14.89 
14.96 
15.04 


>14.00 


13.97 
14.09 
14.18 


13.09 
13.17 
13.25 


13.52 
13.56 
13.61 


< 13.34 


11.76 
12.25 
12.47 


WD 1017-138 


>14.10 


12.93 
13.10 
13.23 


13.05 
13.14 
13.21 


15.73 
15.95 
16.40 


14.46 
14.55 
14.64 


>13.90 


13.90 
14.06 
14.17 


13.02 
13.17 
13.28 


13.55 
13.60 
13.65 




12.18 
12.23 
12.57 


WD 1202+608a 


>14.83 


12.28 
12.60 
12.78 


>14.14 


14.57 
14.70 
14.84 


13.12 
13.37 
13.42 


>14.38 


14.05 
14.11 
14.17 


11.88 
12.18 
12.55 


13.68 
13.72 
13.77 






WD 1202+608b 


>14.65 


12.94 
13.07 
13.17 


>14.14 


15.63 
15.65 
15.68 


14.75 
14.77 
14.79 


>14.10 


14.13 
14.19 
14.25 


12.90 
13.07 
13.19 


13.70 
13.74 
13.78 






WD 1211+332 


>14.39 


13.07 
13.14 
13.19 


>13.49 


15.69 
15.74 
15.79 


14.73 
14.77 
14.81 


>14.03 


14.32 
14.36 
14.40 


13.13 
13.16 
13.19 


13.50 
13.54 
13.58 


<13.43 


12.50 
12.57 
12.63 


WD 1234+481 


>14.23 


12.74 
12.89 
13.01 


>13.60 


15.56 
15.63 
15.71 


14.58 
14.62 
14.67 


>14.15 


13.65 
13.87 
14.01 


12.83 
12.90 
12.96 


12.92 
13.07 
13.19 




12.14 
12.48 
12.66 


WD 1254+223 


13.62 
13.64 
13.66 


<12.66 


12.20: 


14.26 
14.38 
14.58 


13.25 
13.30 
13.45 


12.96 
13.07 
13.16 


<12.69 


11.71 

12.04 
12.23 


<12.54 


<12.94 


<12.30 


WD 1314+293 


14.66 
14.83 
14.92 


<12.73 


12.54 
12.59 
12.64 


14.48 
14.51 
14.54 


13.48 
13.51 
13.54 


13.54 
13.62 
13.69 


12.74 
12.85 
12.94 


11.71 
11.96 
12.11 


12.15 
12.17 
12.19 


<12.64 


<11.97 


WD 1528+487a 


13.26 
13.32 


<12.88 


>13.55: 


<14.07 


<13.46 


<13.22 


<13.04 


<12.55 


<13.22 


<13.30 


<12.56 



13.36 







12 


.89 




15 


.72 


14 


95 




13.80 


13 


.42 


13 


.24 




12 


.10 


WD 1528+487b 


>14.23 


12 


.99 


>13.55: 


15 


.80 


14 


.99 


>14.03 


13.91 


13 


.46 


13 


.36 


<13.30 


12 


.40 






13 


.05 




15 


.88 


15. 


04 




14.01 


13 


.49 


13 


.45 




12 


.70 






13 


.08 


12.78 


15 


.68 


14 


93 


13.61 




13 


.27 


13 


.44 




12 


.39 


WD 1615-154 


>14.03 


13 


.17 


12.96 


15 


.78 


15 


.03 


13.71 


<13.94 


13 


.29 


13 


.50 


<13.16 


12 


.45 






13 


.25 


13.09 


15 


.89 


15 


.18 


13.79 




13 


.31 


13 


.57 




12 


.59 






13 


.02 




15 


.81 


15 


.18 




14.01 


13 


.45 


13 


.20 




12 


.59 


WD 1631+781 


>14.32 


13 


.17 


>13.41 


15 


.90 


15 


28 


>14.12 


14.11 


13 


.53 


13 


.23 


<12.68 


12 


.71 






13 


.32 




15 


.99 


15. 


.38 




14.21 


13 


.61 


13 


.26 




12 


.83 




15.65 








15 


.48 


14 


.58 


14.72 


13.59 


13 


.20 


13 


.01 




12 


.01 
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Table 4 — Continued 



Star 


ClI 


ClI* 


cm 


Ol 


Nl 


Nil 


Sill 


Ar I 


Foil 


FeUI 


PlI 


WD 1634-573 


15.90 
16.25 




<15.4 


15.51 
15.54 


14.62 
14.66 


14.90 
15.08 


13.76 
13.89 


13.26 
13.32 


13.08 
13.15 




12.08 
12.15 


WD 1636+351 


>14.25 


13.23 
13.34 
13.44 


>13.27 


15.58 
15.71 
15.95 


> 15.08 


>14.08 


14.15 
14.30 
14.42 


13.64 
13.70 
13.76 


13.94 
13.98 
14.02 


<13.41 


12.70 
12.75 
12.80 


WD 1800+685 


>14.47 


13.57 
13.63 
13.67 


>13.56 


16.02 
16.12 
16.26 


15.02 
15.05 
15.09 


>14.21 


14.44 
14.53 
14.60 


13.81 

13.83 
13.93 


13.87 
13.97 
14.07 


12.14 
12.80 
13.05 


12.80 
12.94 
13.05 


WD 1844-223 


>14.16 


12.87 
13.04 
13.16 


>13.25 


15.85 
15.97 
16.15 


14.91 
15.06 
15.36 


>13.82 


14.28 
14.38 
14.45 


13.54 
13.58 
13.62 


13.42 
13.53 
13.66 


<13.11 


12.45 
12.59 
12.69 


WD 2004-605 


>14.02 


12.81 
12.99 
13.17 


13.12 
13.15 
13.17 


15.57 
15.65 
15.73 


14.76 
14.84 
14.92 


>13.85 


13.92 
14.07 
14.22 


13.18 
13.26 
13.34 


13.18 
13.25 
13.32 


<12.90 


12.30 
12.37 
12.44 


WD 2011+395 


>14.39 


13.35 
13.38 
13.41 


<13.48 


16.01 
16.06 
16.12 


15.20 
15.48 
15.61 


>14.34 


14.16 
14.23 
14.29 


13.40 
13.43 
13.93 


13.69 
13.72 
13.75 


<14.09 


12.39 
12.59 
12.72 


WD 2111+498 


>14.12 


13.05 
13.12 
13.18 


<13.00 


15.00 
15.06 
15.11 


13.97 
14.00 
14.08 


> 13.91 


13.31 
13.35 
13.31 


12.46 
12.58 
12.68 


13.29 
13.23 
13.35 




<12.47 


WD 2124-224 


> 14.57 


13.41 
13.47 
13.52 


>13.75 


15.91 
15.94 
15.97 


14.93 
14.97 
15.04 


>14.32 


14.44 
14.49 
14.53 


13.50 
13.57 
13.64 


14.18 
14.22 
14.27 




12.62 
12.74 
12.83 


WD 2148+286 


>14.43 


13.73 
13.75 
13.76 


>13.54 


16.17 
16.22 
16.27 


15.48 
15.55 
15.62 


>14.13 




>13.60 








WD 2152-548 


>14.30 


12.65 
12.95 
13.13 


>13.37 


15.42 
15.49 
15.59 


14.37 
14.42 
14.48 


>14.12 


13.38 
13.76 
13.97 


12.56 
12.72 
12.90 


13.34 
13.37 
13.40 


<13.23 


<12.84 


WD 2211-495 


> 14.34 


12.79 
12.88 
12.96 


>13.58 


15.31 
15.34 
15.38 


14.27 
14.30 
14.33 


14.45 
14.60 
14.75 


13.92 
14.00 
14.08 


12.77 
12.82 
12.87 


13.17 
13.25 
13.33 




11.95 
12.05 
12.15 


WD 2247+583 


> 14.41 


13.56 
13.59 
13.63 


> 13.31 


16.53 
16.72 
16.87 


15.66 
15.73 
15.80 


>13.95 


> 14.52 


> 13.50 


13.96 
14.03 
14.10 


<13.14 




WD 2309+105 


>14.20 


13.05 
13.08 
13.10 


>13.37 


15.63 
15.67 
15.71 


14.72 
14.75 
14.78 


>13.81 


14.05 
14.07 
14.10 


13.09 
13.14 
13.21 


13.25 
13.30 
13.35 




12.24 
12.29 
12.34 


WD 2331-475 


>14.43 


13.18 
13.27 
13.36 


>13.60 


15.42 
15.48 
15.54 


14.48 
14.53 
14.58 


>14.27 


13.89 
13.99 
14.09 


12.67 
12.77 
13.87 


13.27 
13.32 
13.37 




12.08 
12.18 
12.28 



Note. — For a given sightline and a given species, there are 3 entries: the second number corresponds to the most likely value, 
while the first and third correspond to — lcr and + 1ct, respectively. < indicates a 3tr upper limit detection, except in the case of 
ClII where the upper limit reflects that this line may be blended with the C III photosphcric lines. > indicates a lower limit, i.e. 
the line is saturated. Colon means that the value is uncertain. Sec § 3 and Table 1 for more details and references, respectively. 
Two entries for one sightline (denoted by a and b) indicate that two components are detected with FUSE: for WD 0455 — 282, 
component a is at —58 kms -1 and b at —2 kms -1 ; for WD 1202+608, component a is at —58 kms -1 and b at —2 kms -1 ; and 
for WD 1528+487, component a is at —56 kms -1 and b at —16 kms -1 . 



Table 5 
h 2 measurements 



Sightline 




= 






AT. 7=1 






Nj =2 




Nj =3 




Wall 




b 




T cx 






























(km s 




(K) 


WD 0004+330 


13.45 


± 


0.06 


(1) 


14.35 


± 0.04 


(22) 


13.71 


± 0.04 


(10) 


13.52 ± 0.07 


(4) 


14.46 


3.8 


± 


0.7 


314 ± 42 a 


WD 1636+351 


14.21 


± 


0.20 


(7) 


14.83 


± 0.07 


(25) 


14.30 


± 0.15 


(13) 


14.04 ± 


(3) 


15.05 


3.1 


± 


0.2 


314 ± 31 


WD 1800+685 


14.02 


± 


0.26 
0.10 


(5) 


14.49 


± 0.15 


(12) 


14.00 


± 0.10 


(12) 


13.79 ± 0.18 


(1) 


14.75 


5.0 


± 


4.2 

u 


301 ± 38 


WD 2011+395 










14.55 


± 0.10 


(13) 


14. OF. 


. ± 0.20 


(2) 


13.83 ± 0.20 


(2) 


14.22 


2.8 


± 


0.3 


349 ± 74 


WD 2148+286 


14.32 


± 


0.11 
0.09 


(4) 


15. 2E 


_i_ 0.24 
' ^ 0.12 


(8) 


14.34 


i 0.15 
1 0.16 


(13) 


14-18 ± g;i; 


(14) 


15.38 


3.0 


± 


0.6 


297 ± 35 


WD 2247+583 


14.24 


± 


0.03 


(6) 


14.94 


± 0.04 


(13) 


14.23 


± 0.03 


(14) 


13.94 ±0.04 


(7) 


15.11 


4.1 


± 


0.4 


276 ± 8 



Note. — The number of lines used to measure the column density is shown in parentheses. For WD 0004+330, WD 2148+286, 
and WD 2247+583, the column densities and 6-values are from Oliveira et al. (2003, private communication) Sonncborn et al. (2003, 
erratum in preparation), and Oliveira et al. (2003), respectively, a : For WD 0004+330, only the J-levels 1,2,3 could be fitted in the 
temperature relation within la error bar, Toi ~ 1363 K. 
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Table 6 

[ArI/O I], [N I/O I], AND N I/(N II + N I) RATIOS 



Sightlines 



[Ar I/O I] 



[N I/O I] N I/(N II+N I) 



WD 0004+330 
WD 0050-332 
WD 0131-163 
WD 0455-282b 
WD 0501+527 
WD 0549+158 
WD 0621-376 
WD 0715-703 
WD 1017-138 
WD 1202+608a 
WD 1202+608b 
WD 1211+332 
WD 1234+481 
WD 1254+223 
WD 1314+293 
WD 1528+487b 
WD 1615-154 
WD 1631+781 
WD 1634-573 
WD 1636+351 
WD 1800+685 
WD 1844-223 
WD 2004-605 
WD 2011+395 
WD 2111+498 
WD 2124-224 
WD 2148+286 
WD 2152-548 
WD 2211-495 
WD 2247+583 
WD 2309+105 
WD 2331-475 



> -0.41 

< -0.45 
-0.78+ 

< -0.15 

< -0.72 
< 0.33 

-0.26 ± H-S| 
-0.59 ± «;| 

-°- 64± III 
-0.38 ± «;| 

— 44 + °- 12 

— 44 + 0.06 
U - 44± Q-05 

-0.59 ± 

-0.20 ± 
-0.41 ± 
-0.20 ± 
-0.35 ± 
-0.23 ± 
-0.11 ± 
0-13 ±8;i| 
-0.15 ± 
-0.25 ± 
-0.25 ± 
-0.48 ± 
-0.34 ± 
-0.23 ± 

> -0.48 
-0.63 ± 
-0.38 ± 
-0.65 ± 
-0.39 ± 
-0.64 ± 



Mi 

15 

m 
a 

IS 
4? 

0„Q7 



8:10 
8:1? 
8:48 

m 

0.04 



He 



-0.29+ 0-20 

-0.20 ± 

-0.22 ± 

-0.71""" 
_n 95 4- 0.06 

u - zc> ^ 0.06 

< 0.09 

-0.16 ±8:8| 

-0.18 ± 
-0.64 + 
-0.57 + 
-0.12 ± 
-0.21 ± 
-0.25 ± 
-0.32 ± 
-0.24 + 
-0.05 ± 

0.01 ± 

0.14 + 

u - ld 33 0.05 
> 0.13 

"°- 31 ± Ui 
-0.15+ ; 

-0.05 ± J];" 

0.19 ± Sis, 
-0.30 ±8;| 

-0.20 + 



8:i§ 

8:?e 
8:81 
8:8§ 

ni 
m 

8:41 

m 

8:11 

0J5 



0.09 ± 8-8 



8:?4 
0.10 



-0.31 ± 
-0.28 ± 
-0.23 ± 
-0.16 ± 
-0.31 ± 



i:8e 



< 
< 
< 
< 
< 
0. 
< 
< 
< 
< 
< 
< 
< 

0.63 

0.44 
< 

0.95 
< 

0.34 
< 
< 
< 
< 
< 
< 
< 
< 
< 

0.33 
< 
< 
< 



0.10 



0.94 
0.71 
0.86 
0.25 
0.52 
56 : 
0.77 
0.90 
0.82 
0.09 
0.82 
0.85 
0.75 

+ »:iin 

33 0.03 

0.90 

4. 0.08 
33 0.04 

0.94 

_l_ 0.13 
33 0.09 

0.91 
0.87 
0.94 
0.91 
0.93 
0.55 
0.82 
0.96 
0.67 

_i_ 0.08 
31 0.06 

0.98 
0.89 
0.63 



Note. — Suffix "a" after star names means it is a high velocity 
component (see § 9), while "b" means that it is the low velocity com- 
ponent. Colon means that the value is uncertain. 



Table 7 



C III/C II RATIO 


Targets 


log[JV(C III/AT(C II)] 


log[JV(C III/iV(C II)] C 


Note 


WD 0050-332 


< -0.95 


< -1.26 to < -0.95 


C II LL, C III SC? 


WD 0131-163 


< -0.93 


< -0.93 


C II LL 


WD 0455-282a 


< -0.92 




C III UL 


WD 0549+158 


-1.09 


< -1.09 to < -1.05 


C III SC? 


WD 0621-376 


< -0.90 


< -1.05 


C II LL, C III SC 


WD 1017-138 


< -0.96 


< -0.96 


C II LL 


WD 1254+223 


-1.44 : 


< -1.44 : 


C III SC? 


WD 1314+293 


-2.24 


-2.39 to -2.24 


C III SC? 


WD 1528+487a 


> 0.23 : 




C III LL 


WD 1615-154 


< -1.07 


< -1.07 


C II LL 


WD 1634-573 


< -0.50 


-0.63 


C III SC 


WD 2004-605 


< -0.87 


< -0.87 


C II LL 


WD 2011+395 


< -0.91 


< -1.44 


C II LL, C III SC 


WD 2111+498 


< -1.12 


< -1.52 


C II LL, C III SC 



Note. — The second column presents the raw measurement, while the third one includes 
possible stellar correction for C III. When two numbers are indicated in the third column, it 
means that only an upper limit was derived for the amount of photospheric C III A977. LL: 
lower limit. UL: upper limit. SC: C III stellar contamination. SC?: possible C III stellar 
contamination. Suffix "a" after star names means the ratio is for a high velocity component 
(see § 9), and this component is not blended with any stellar lines. Colon means that the 
value is uncertain. 
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Table 8 
Fe III/Fe II RATIO 



Targets 


log[7V(Fe III/7V(Fc II)] 


WD 0004+330 


< -0.65 


WD 0455-282b 


< 0.05 


WD 0715-703 


< -0.22 


WD 1211+332 


< -0.11 


WD 1314+293 


< 0.47 


WD 1528+487b 


< -0.06 


WD 1615-154 


< -0.34 


WD 1631+781 


< -0.63 


WD 1636+351 


< -0.57 


WD 1800+685 


-1.17 : 


WD 1844-223 


< -0.42 


WD 2004-605 


< -0.35 


WD 2011+395 


< 0.37 


WD 2152-548 


< -0.14 


WD 2247+583 


< -0.89 



Note. — Here "<" is a 3cr upper limit. 
Colon means that the value is uncertain. 
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Fig. 1. — Target locations and the total O I column densities in Galactic coordinates (l,b). The diameter of each circle is inversely 
proportional to the distance of the line of sight, and the shading of the symbol indicates the total column density. The H I values were 
computed from the O I values using the O I/H I ratio discussed in the text. We indicate the presence of molecular hydrogen by the symbol 
H2 for detection by this survey and "C" for sightlines near WD 1615— 154 and I ~ 150° from the survey of Spitzer, Cochran, & Hirshfeld 
(1974) based on Copernicus data (the stars in the Copernicus study are between 111 and 166 pc). 
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Fig. 2. — Two examples of FUSE WD spectra, with WD 0131—163 shown in top three panels and WD 1800+685 in bottom three panels. 
WD 1800+685 is a pure hydrogen WD with no photosphcric metal absorption lines, but shows absorption by interstellar atoms, ions, and H2. 
WD 0131— 163 has a few stellar lines (note in particular the stellar Si IV blended with Ar I at 1066.66 A). The positions of strong airglow 
emission lines are indicated by the earth symbol. 
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Fig. 3. — Histograms of column densities for the different ions observed. The mean, standard deviation and number of detected clouds (n) 
are listed in the upper-left of each plot. 
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Fig. 4. — Logarithmic ratios of Ar I, N I, and (N I + N II)/0 I plotted versus the total O I column density and the distance of the sightlincs 
(note that the distance is only an upper limit to the interstellar clouds). The horizontal lines indicate the solar, B-type star, and interstellar 
values. The Local Bubble (log N (O I) < 15.60) data points are shown with filled circles. 
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Fig. 5. — Logarithmic ratio of Ar I/O I plotted versus N I/O I. Both ratios arc normalized to their solar values. The straight continuous 
line shows a 1:1 relationship. The dashed lines present the expected deficiency of Ar I and N I from a photoionization model (see § 6.3 and 
Figure 6). The Local Bubble (log7V(0 I) < 15.60) data points are shown with filled circles. 
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Fig. 6. — Results from a photoionization model of gas in the LB. Top: Logarithms of the expected deficiencies in the neutral fractions of 
fo = ra(ncutral)/n(total) of N, 0, and Ar compared to those of H (see 6.3 for more details). Also plotted are the logarithms of n(X+)/n(total), 
i.c the singly ionized species, compared to the neutral fraction of H, along with the ratios C III/C II and Fe III/Fe II (using the same 
logarithmic scale on the vertical axes) . All the calculations assume the same ionizing radiation field and include the effects of charge exchange 
with neutral hydrogen These fractions are plotted as a function of absorbing column and hence shielding depth scaled to N(H I) for a cloud 
of uniform pressure p/k. Bottom: Densities of neutral hydrogen, neutral helium, and electrons against shielding depth. 
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Fig. 7. — Logarithmic ratios of of Fe II and P II to Si II against the total O I column density and against the distance of the 
The dashed lines indicate the solar value of those ratios. The Local Bubble (log7V(0 I) < 15.60) data points are shown with filled 



sightlines. 
circles. 
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Fig. 8. — Logarithmic ratio of P Il/Si II against Fc II/P II normalized to their solar values. The "Cold", "Warm", and "Halo" refer to 
depletion patterns observed in our Galaxy. We also indicate for each point the first 4 digits of the WD. 
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Fig. 9. — Logarithmic of the cooling rate in ergs -1 per H I atom (upper panel) and of the electron density in cm~ 3 (lower panel) against 
the total H I column density (the latter being derived from O I: log ^(H I) = log iV(0 I) + 3.5) (see § 8 for more details). For clarity, we do 
not show the lower limits, but none contradicts the measured values. Also shown as a solid line is logn e from Figure 6. 



